Mon. Not. R. Astron. Soc. 000, [T]-?? (2008) Printed 25 September 2008 



(MN JATeH style file v2.2) 



Galaxy Merger Morphologies and Time-Scales 

from Simulations of Equal-Mass Gas-Rich Disc Mergers 



oo 

o 
o 

(N 

in 

(N 

6 

ti 



(N 
> 

(N 

in 
o 
oo 
O 



X 



Jennifer M. Lotz/*, Patrik Jonsson,^ T.J. Cox,^f and Joel R. Primack^ 

^ National Optical Astronomical Observatory, 950 N. Cherry Avenue, Tucson, AZ 85719, USA; lotz@noao.edu 

^ Department of Physics, University of California, Santa Cruz, OA, USA, 95064 

^ Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA, 02138 USA 



Submitted 2008 8 May 



ABSTRACT 

A key obstacle to understanding the galaxy merger rate and its role in galaxy evo- 
lution is the difficulty in constraining the merger properties and time-scales from in- 
stantaneous snapshots of the real universe. The most common way to identify galaxy 
mergers is by morphology, yet current theoretical calculations of the time-scales for 
galaxy disturbances are quite crude. We present a morphological analysis of a large 
suite of GADGET N-Body/hydro-dynamical equal-mass gas-rich disc galaxy mergers 
which have been processed through the Monte-Carlo radiative transfer code SUNRISE. 
With the resulting images, we examine the dependence of quantitative morphology (G, 
M20, C, A) in the SDSS g-band on merger stage, dust, viewing angle, orbital parame- 
ters, gas properties, supernova feedback, and total mass. We find that mergers appear 
most disturbed in G — M20 and asymmetry at the first pass and at the final coales- 
cence of their nuclei, but can have normal quantitative morphologies at other merger 
stages. The merger observability time-scales depend on the method used to identify 
the merger as well as the gas fraction, pericentric distance, and relative orientation 
of the merging galaxies. Enhanced star formation peaks after and lasts significantly 
longer than strong morphological disturbances. Despite their massive bulges, the ma- 
jority of merger remnants appear disc-like and dusty in g-band light because of the 
presence of a low-mass star-forming disc. 
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1 INTRODUCTION 

It remains unknown to what degree present-day galaxies are 
assembled discretely via the merger of pre-existing galaxies 
(e.g. KaufFmann, White, & Guideroni 1993; Somerville, Pri- 
mack, & Faber 2001) or through more continuous processes 
such as cold gas and dark matter accretion (e.g. Keres et 
al. 2005; Dekel & Birnboim 2006). The most obvious way 
to constrain the importance of galaxy mergers is to count 
the number of on-going merger events. However, current ob- 
servational constraints on the galaxy merger rate are highly 
uncertain and strongly debated (Ryan et al. 2008; Lin et al. 
2008; Lotz et al. 2008; Renzini 2007; Kartaltepe et al. 2007; 
Masjedi et al. 2008). Moreover, theoretical predictions for 
the galaxy merger rate and mass assembly can vary by fac- 
tors of ten (Jogee et al. 2008). These discrepancies are par- 
tially the result of the non-trivial conversion of the observed 
number density of galaxy mergers into a galaxy merger rate 
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and the comparison of this galaxy merger rate to the cos- 
mological predictions for dark matter halo assembly (e.g. 
Berrier et al. 2007; Guo & White 2008; Kitzbichler & White 
2008). 

Since the realization that the merger of two disc 
galaxies could produce a spheroidal galaxy (Toomre 1977), 
spheroidal galaxies, red galaxies, and post-starburst galaxies 
have been used to indirectly trace the role of galaxy merg- 
ers in galaxy evolution (e.g. Bell et al. 2004; Faber et al. 
2007; Brown et al. 2007; Abraham et al. 2007; Hopkins et 
al. 2007; Hogg et al. 2006). However, there are multiple ways 
to quench star-formation and produce post-starbursts or red 
spheroidal galaxies (e.g. Dekel & Birnboim 2006; Moore, 
Lake, & Katz et al. 1998). Also, simulations of galaxy merg- 
ers often show significant star-formation and disc compo- 
nents well after the merger event (Springel & Hernquist 
2005; Robertson et al. 2006a; Naab, Jesseit, & Burkert 2006; 
Cox et al. 2006, 2008; Khalatyan et al 2008), thus the photo- 
metric and morphological signatures of merger remnants are 
ambiguous. A merger remnant's star-formation history, mor- 
phology, and kinematics are likely to depend on the proper- 
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ties of the progenitors and the merger conditions, which arc 
increasingly difflcuh to determine as time passes. Therefore 
direct observations of on-going galaxy mergers are needed to 
constrain the role of mergers iu galaxy assembly. 

Galaxy merging is a process that lasts several billion 
years rather than a short-lived event. The signatures of a 
galaxy merger change with merger stage, making the identi- 
fication of galaxy mergers challenging. If hierarchical models 
of galaxy assembly are correct, then the majority of massive 
galaxies could be considered an on-going merger or a merger 
remnant which will undergo another merger event within 
a few Gyr (Stewart et al. 2008). We shall define a galaxy 
merger as a pair of galaxies which arc gravitationally bound 
and whose orbits will dynamically decay such that their nu- 
clei will merge within x billion years, where x is typically 
1-3 Gyr for 'major' mergers with mass ratios greater than 
1:3. 

There arc two general approaches to identifying such 
systems observationally. The first approach is to find close 
pairs of galaxies before their nuclei have coalesced, either by 
selecting galaxies close in projected angular separation and 
line-of-sight radial velocity (e.g. Patton et al. 2000; Barton, 
Geller, & Kcnyon 2000; Lin et al. 2004; de Propris et al. 
2005; Lin et al. 2008) or by measuring the deprojected cor- 
relation function on small scales (Bell et al. 2006a; Masjedi 
et al. 2006, 2008; Li et al. 2008). The second approach is to 
identify morphologically disturbed galaxies, some of which 
will bo post-mergers and some of which will be interacting 
pairs. Morphological disturbances can be found qualitatively 
through visual inspection (e.g. Kampczyk et al. 2007; Bundy 
et al. 2005; Brinchmann et al. 1998), or by quantitative mea- 
sures such as the Gini coefficient, second-order moment of 
the brightest 20% of the light (M20), and asymmetry (Abra- 
ham et al. 1994; Conselice 2003; Abraham et al. 2003; Lotz, 
Primack, & Madau 2004 [LPM04]; Scarlata et al. 2007). At 
present, all quantitative merger indicators are calibrated em- 
pirically using galaxies with 'normal' and 'disturbed' visual 
classifications. 

Translating the number of observed merger candidates 
into a merger rate requires the assumption of an observ- 
ability time-scale - the time during which one would have 
identified the system as merging. Until now, this time-scale 
has been poorly constrained. Close pair studies often as- 
sume a dynamical friction time-scale (which varies from 200 
Myr to 1 Gyr). However this value does not take into ac- 
count the range of possible orbits for the merging system 
nor the time during which the system would not meet the 
pair criteria at very large and very small separations. Also, 
recent comparison of N-body simulations to analytical calcu- 
lations indicate that analytically-derived dynamical friction 
time-scales can deviate substantially from those predicted by 
N-body simulations (Boylan-Kolchin, Ma, & Quataert 2008; 
Jiang et al. 2008) . Very similar time-scales are also generally 
assumed for mergers selected using both visual and quanti- 
tative morphologies. This assumption is even less likely to 
be valid given that different morphological selection criteria 
are sensitive to different stages of the merger process. For 
example, visual classification using a combination of signa- 
tures (e.g. tidal tails, multiple nuclei, shells) is likely to be 
sensitive for longer time-scales and lower-mass merger ratios 
than current quantitative methods. 

The observability time-scale for a particular merger may 



depend on (1) the method used to identify the merger; 
(2) the merger parameters (mass ratio, gas properties, 
bulge/disc ratio, orbits, dust content); and (3) the observa- 
tional selection (observed wavelength, viewing angle, spatial 
resolution). Cosmological-scale numerical simulations cur- 
rently do not have the spatial resolution to directly de- 
termine the cosmologically-averaged observability time-scale 
for each method. Therefore, one is required to use a suite of 
galaxy-scale numerical simulations which span a largo range 
of input merger parameters to constrain the observability 
time-scales for the different input parameters. Given a suffi- 
ciently broad range of merger parameters, the observability 
time-scales for each parameter set may then be weighted 
by the probability distribution of the mass ratios, gas frac- 
tion, etc. which can be computed from current cosmological- 
scale simulations. An additional complication is that galaxy- 
scale numerical simulations typically track the distribution 
of 'particles' (star, gas, and dark matter), as opposed to 
the projected light distribution at a particular wavelength 
(which is what is observed) . A few works have attempted to 
quantify the observability time-scales using only the stellar 
particles (Conselice 2006; Bell et al. 2006b) or gas particles 
(lono, Yun, & Mihos 2004) and ignoring the effects of dust 
and age-dependent stellar luminosities. While this may be 
acceptable for gas-poor dissipationless mergers (e.g. Bell et 
al. 2006b), the appearance of most gas-rich mergers is almost 
certainly affected strongly by both young stars and dust. 

In this work, we present a first attempt to constrain the 
observability time-scales for a variety of methods for iden- 
tifying galaxy mergers. We present a morphological analy- 
sis of a large suite of GADGET N-Body/SPH equal-mass 
gas-rich disc galaxy merger simulations which have been 
processed through the Monte-Carlo radiative transfer code 
SUNRISE. With the resulting images, we examine the de- 
pendence of quantitative morphology in the SDSS 5-band 
(A, = 4686A) on merger stage, dust, viewing angle, orbital 
parameters, gas properties, supernova feedback, and total 
mass. We constrain the time-scales of quantitative morphol- 
ogy disturbances in Gini coefficient, M20, and asymmetry, 
and the time-scales during which close pairs lie at projected 
separations Rproj <20, 30, 50, and 100 kpc, assuming 
h — 0.7. Finally, we compare the simulated merger remnant 
morphologies and star-formation rates. In §2, we describe 
the galajcy merger simulations, including the GADGET N- 
Body/SPH calculations, the SUNRISE Monte-Carlo radia- 
tive transfer calculations, the initial galaxy models, and the 
range of merger parameters explored. In §3, we define the 
morphological quantities G, M20, C and A. We also describe 
the different observational criteria used to identify galax:y 
mergers, and define the merger observability time-scale for 
each method. In §4, we present the results of our analysis, 
and in §5 we discuss the implications of these results for 
finding galaxy mergers, calculating the merger rate, and the 
properties of merger remnants. A subsequent paper with a 
similar analysis of non-equal mass mergers is in preparation. 

2 GALAXY MERGER SIMULATIONS 

2.1 GADGET N-Body/SPH simulations 

AH of the simulations presented in this work were per- 
formed using the N-Body/SPH code GADGET (Springel 
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et al. 2001). The details of these simulations, their global 
star-formation histories, and their remnant properties are 
discussed in Cox et. al (2004, 2006, 2008). Each galaxy is 
initially modeled as a disc of stars and gas, a stellar bulge, 
and a dark matter halo. The stellar and dark matter par- 
ticles are coUisionless and are subject to only gravitational 
forces. The gas particles are also subject to hydro-dynamical 
forces. The baryonic and dark matter particles have gravi- 
tational softening lengths of 100 pc and 400 pc respectively. 
The SPH smoothing length for the gas particles indicates 
the size of the region over which the particle's hydrodynamic 
quantities are averaged and is required to be greater than 
half the gravitational softening length or > 50 pc. While we 
use the first version of GADGET (Springcl et al. 2001), the 
smoothed particle hydrodynamics are upgraded to use the 
'conservative entropy' version that is described in Springel 
& Hernquist (2002). The radiative cooling rate A„et(p, u) is 
computed for a primordial plasma as described in Katz et 
al. (1996). 

Gas particles are transformed into coUisionless star par- 
ticles assuming the Kennicutt-Schmidt law (Kennicutt 1998) 
where the star-formation rate depends on the local gas den- 
sity Pgas- This occurs in a stochastic sense (Springel & Hern- 
quist 2003) in which each gas particle can spawn one or two 
new star particles with a probability determined by the star- 
formation rate. These new star particles have typical masses 
~ 10'' Mq, and are assigned ages based on their formation 
time and metallicities based on the metallicity of the gas 
particle from which they are spawned. We adopt the in- 
stantaneous recycling approximation for metal production 
whereby massive stars arc assumed to instantly become su- 
pernovae, and the metals produced are put back into the 
gas phase of the particle. In this version of GADGET, met- 
als do not mix and remain in the gas particle in which they 
are formed. The enriched gas contribution from stellar winds 
and Type la supernovae are ignored. Unlike the metals, there 
is no recycling of hydrogen and helium to the gas. 

Feedback from supernovae is required to produce sta- 
ble star-forming discs. Energy from supernovae heats and 
pressurizes the interstellar gas and stabilizes it against grav- 
itational collapse. Because of the limited resolution of most 
N-Body/SPH simulations, the physical processes associated 
with supernova feedback cannot be directly modeled and 
must be included using simple prescriptions. We test two 
supernova feedback models, which are discussed in detail in 
Cox et al. (2006). Both models assume the supernova feed- 
back energy is dissipated on an 8 Myr time-scale, and have 
a equation of state parameterized by n, where P ~ p^ts"^^^ ■ 
The n = 2 model treats star-forming gas with a stiff equa- 
tion of state where the pressure in star-forming regions scales 
as P ~ pgas, while the n = model assumes that this 
gas is isothermal with an effective temperature ~ 10® K 
{P Pgas)- Both feedback models produce stable isolated 
star-forming discs and predict similar gas consumption dur- 
ing the merger, but the strength and timing of the merger- 
induced starbursts depend on the feedback model assumed 
(Cox ct al. 2006). No active galactic nuclei arc included 
in these simulations. Such AGN are expected to influence 
the star-formation and morphologies only during the post- 
merger stages (see §5 for discussion). 



2.2 SUNRISE Monte Carlo radiative transfer 
processing 

SUNRISE is a parallel code which performs full Monte Carlo 
radiative-transfer calculations using an adaptive-mesh re- 
finement grid (Jonsson 2006; Jonsson et al. 2006). This code 

was developed to calculate the effects of dust on the emis- 
sion from the GADGET N-Body/SPH simulations. SUN- 
RISE can model arbitrary geometries of emitting and ab- 
sorbing/scattering material with a large spatial dynamical 
range and efficiently generate images of the emerging radia- 
tion viewed from arbitrary points in space. We use SUNRISE 
v2 for these simulations (Jonsson 2006). Given a particular 
simulation geometry and viewing angle, SUNRISE v2 per- 
forms the Monte-Carlo radiative transfer calculation for 20 
wavelengths from the far-ultraviolet to the mid-infrared and 
interpolates a resulting spectral energy distribution of 510 
wavelengths including the effects of absorption and scatter- 
ing. 

At least 30 timesteps are analysed for each merger sim- 
ulation. For each GADGET simulation timestep, SUNRISE 
assigns a spectral energy distribution to each star particle 
using the STARBURST99 population synthesis models (Lei- 
therer et al. 1999). New star particles are treated as single 
stellar populations with ages based on their formation time. 
Bulge star particles form in an instantaneous burst well be- 
fore the start of the simulation (8-13 Gyr depending on the 
galaxy model; see Rocha et al. 2008). Initial disc star parti- 
cles are assumed to have formed with an exponentially de- 
clining star-formation starting at the time of the formation 
of the bulge. The metallicities of the gas and stars of the 
initial galaxy models decline exponentially with the radius 
of the disc. The density of dust is linearly proportional to 
the density of metals in the gas. 

The dust attenuation in the initial galax;y models have 
been found to reproduce the observed galaxy inclination- 
attenuation relations and the global infrared-to-ultraviolet 
flux ratios for spiral galaxies (Rocha et al. 2008). During 
the merger, the gas discs of the initial galaxies become dis- 
rupted, resulting in complicated dust geometries. The full 
radiative transfer calculation done by the SUNRISE v2 code 
are well-suited to determining the effects of complicated 
dust geometries. However, the dust approximation is lim- 
ited by the low spatial resolution of the regions of cold gas 
and star-formation in the input GADGET simulations (~ 
100 pc). Therefore the atteimation of very young stars may 
be underestimated, and the effects of clumpy dust and gas 
on small scales are not included. A future version of SUN- 
RISE will include improved treatment of small scales. The 
attenuation and infrared luminosities of the output SUN- 
RISE images during the merger have been compared to the 
available literature for local dusty mergers. Jonsson et al. 
(2006) found that these simulations reproduce the observed 
relationships between ultraviolet spectral slope /3 and the 
global infrared-to-ultraviolet flux ratio, thus the simulated 
dust distribution is a reasonable approximation of the dust 
found in local dusty mergers. 

Images in multiple band-passes (GALEX FUV/NUV, 
SDSS ugriz, 2MASS JHK) for 11 isotropically positioned 
viewpoints ('cameras') are generated and the total absorbed 
bolometric luminosity over all wavelengths/viewing angles is 
calculated. These cameras are positioned with respect to the 
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Figure 1. SDSS u — r — z composite colour images with dust extinction for the high-resolution Sbc prograde-prograde simulation 
(SbcPPxlO) as viewed by camera (face-on). Time since the start of the simulation in given in the upper left corner of each image. 
Shown in the top row are the initial pre-merger galaxies, the first pass, the maximal separation after the first pass, and in the bottom 
row are the merger of the nuclei, the post-merger at 0.5 Gyr after the merger, and the remnant at 1 Gyr after the merger. The field of 
view for the initial galaxies and the maximal separation is 200 kpc, while the field of view for the other images is 100 kpc. The merger 
morphologies are most disturbed at the first pass and merger. Star-forming regions in the initial discs, tidal tails, and outer regions of 
the remnant appear blue, while the dust-enshrouded star-forming nuclei appear red. 



plane of the merger orbit at (j), 6 (degrees) = (0,0), (0, 79), 
(72, 79), (144, 79), (216,79), (288,79), (0, 127), (72, 127), 
(144, 127), (216, 127), and (288,127). In Figures 1 and 2, we 
show examples of composite SDSS u — r — z images for one 
of our simulations viewed face-on (camera 0) and edge-on 
(camera 4). The predicted dust attenuation for the initial 
undisturbed galaxy models agrees well with observations of 
dust attenuation in local disc galaxies (Rocha et al. 2008). 

2.3 Initial Galaxy Models 

The goal of this work is to calibrate the morphological dis- 
turbance time-scales for merging and interacting galaxies 
using realistic 'observations' of galaxy merger simulations in- 
cluding the effects of star formation and dust. Dust and star 
formation have a much stronger effect on gas-rich merger 
morphologies than dissipationless mergers, hence we exam- 
ine the mergers of gas-rich disc galaxies. All of the disc 
galaxy models explored here have relatively small bulge com- 
ponents with stellar bulge-to-disc mass ratios < 0.25. Such 
low bulge-to-disc ratio galaxies are more likely to experience 
strong starbursts (Mihos & Hernquist 1996, Cox et al. 2008) 
and stronger morphological disturbances (Conselice 2006). 
These galaxies may also be more representative of high red- 
shift mergers, as bulge-dominated systems are increasingly 



rare at 2; > 1 (e.g. Lotz et al. 2008; Ravindranath et al. 
2006). The structure of dissipationless spheroidal galaxy 
merger remnants have been studied by several other authors 
(e.g. Naab, Khochfar,& Burkert 2006; Boylan-Kolchan, Ma, 
& Quartaert 2005). 

We adopt two general models for gas-rich discs: the 
'Sbc' model tuned to match a large, gas-rich Sbc disc galaxy, 
and the 'G-series' of discs with varying mass (G3, G3, Gl, 
GO) with lower gas fractions tuned to match SDSS obser- 
vations of local galaxies. (We use the notation of 'Sbc' and 
'G' for these different galaxy models in keeping with pre- 
vious publications based on these simulations: Cox et al. 
2006, 2008, Jonsson et al. 2006, Rocha et al. 2008). Each 
galaxy model contains a rotationally supported disc of gas 
and stars, a non-rotating stellar bulge, and a massive dark 
matter halo (Table 1). A detailed description of the galaxy 
disc models can be found in Cox et al. (2006, 2008), Jonsson 
et al. 2006 and Rocha et al. (2008). 

The Sbc model parameters are motivated by observa- 
tions of local gas-rich, disc-dominated Sbc galaxies similar 
to the Milky Way (see Table 1; Cox et al. 2006). The optical 
(stellar) disc scalelength, dynamical mass, and gas fraction 
are taken from Roberts & Haynes (1994). The bulge-to-disc 
ratio and bulge size are from the observations of de Jong 
(1996) and the total stellar mass is derived using the Bell 
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Figure 2. SDSS u — r — z composite colour images for tlie same simulation as Fig. 1 (SbcPPxlO) as viewed by camera 4 (roughly 
edge-on). Time since the start of the simulation in given in the upper left corner of each image. The timesteps and image scales are same 
as the previous figure. When viewed edge-on, the dust lanes associated with initial discs and remnants are clearly visible. 



Table 1. Initial Galaxy Conditions 
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2.0- 10" 


12 


7.0- 10^ 


4.7 • 10^ 


3.0 ■ 10** 


2.0 • 10" 


0.04 


0.29 


1.48 


0.33 
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0.01 
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Total number of particles in GADGET simulation for fiducial resolution models. 
Virial mass 

Dark matter halo concentration 
Baryonic mass 
Mass of stellar disc 
Mass of stellar bulge 
Mass of gaseous disc 
Fraction of baryons in the bulge 
Fraction of baryons in gas 
Scalelength of stellar disc 
Scalelength of bulge 
Scalelength of gaseous disc 
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& de Jong (2001) relations. We assume that the gas disc is 
exponential with a scalelength three times the stellar disc 
scalelength (Broeils & van Woerden 1994). We adopt an 
adiabatically contracted NFW dark matter halo with a con- 
centration of 11. The resulting Sbc model has a viral mass of 
8.12 xlO" M0, with a 12.5% baryonic mass fraction. Fifty- 
two percent of the baryons are in gas (mostly at large radii), 
and ~ 10% of the baryons are bulge stars. 

In order to sample the parameter space spanned by 
many present-day galaxies, we also explored mergers be- 
tween model galaxies with masses, bulge-to-disc ratios, and 
gas fractions motivated by SDSS estimates of typical local 
galaxies (Table 1; Cox et al. 2008). We refer to these model 
galaxies as the G-series. The largest galaxy (G3) is chosen 
to have a stellar mass ~ 5 x 10^" Mq , and the smaller galax- 
ies are chosen to have stellar masses ^ 1.5 x 10'" Mo (G2), 
0.5 X 10'° M0 (Gl), and 0.1 x 10^° Mq (GO), spanning a fac- 
tor of 50 in stellar mass. The stellar half-light radii are from 
the stellar mass-size relation of Shen et al. (2003) . The bulge- 
to-disc ratios are taken from de Jong (1996) and used to de- 
termine the stellar disc and bulge masses and scalelengths. 
The gas fractions and masses are determined from the gas 
mass - stellar mass scaling relation from Bell et al. (2003). 
As for the Sbc model, the gas scalelength is assumed to be 
three times the stellar disc scalelength. We adopt NFW dark 
matter halo profiles selected such that the rotation curves lie 
on the baryonic TuUy-Fisher relation (Bell & de Jong 2001; 
Bell et al. 2003) . Unlike the Sbc model, these models do not 
include adiabatic contraction. The total mass-to-light ratio 
is assumed to vary with mass such that lower mass galaxy 
model have higher mass-to-light ratios and the range in to- 
tal mass is a factor of 23. While the total masses of the Sbc 
and G3 models are similar, the G3 model has a lower gas 
fraction, a smaller disc scalelength, and consequently, much 
less gas at large radii than the Sbc model (Table 1). 

2.4 Galaxy Merger Parameters 

Each of the galaxy merger simulations presented here are 
mergers of identical equal-mass galaxies; mergers of non- 
equal mass galaxies will be presented in a later paper. 
The Sbc-Sbc merger simulation parameters are selected 
to probe a range of different merger orientations and or- 
bits (Table 2), including parabolic orbits with roughly 
prograde-prograde (SbcPP), prograde-retrograde (SbcPR), 
retrograde-retrograde (SbcRR), and prograde-polar orien- 
tations (SbcPol), parabolic prograde-prograde mergers with 
very small/large pericentric distances (SbcPPr-, SbcPPr-f), 
and a highly radial orbit with prograde-retrograde orienta- 
tion (SbcR). For most of the simulations, the galaxy orbits 
are initialized to be parabolic with the given pericentric dis- 
tance; as the simulations progress, the galaxies lose angu- 
lar momentum because of dynamical friction and eventually 
merge. The galaxies in the radial orbit simulation (SbcR) 
start out nearly at rest, so the pericentric distance and ec- 
centricity have little physical meaning in this case. The G-G 
merger simulations were all run with the same orientation 
(prograde-prograde) and slightly sub-parabolic orbits (ec- 
centricity = 0.95) to probe the effect of varying total mass 
(G3PP, G2PP, GIPP, GOPP; Table 2). Sub-parabolic orbits 
were chosen so that the lower mass mergers (Gl, GO) would 
merge in less than 2 Gyr. This choice will affect the close 
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" The numerical resolution of simulation is xNpart, where Npart 
is given in Tabic 1. 

^ Supernova feedback polytropic index n where P oc pi+("/2) in 
star-forming regions. 

Initial orientation of galaxy 1 with respect to the plane of the 

orbit in spherical coordinates, where d = arctan( ) and 

(f) = arctan(^) 

Initial orientation of galaxy 2 
^ Eccentricity of the orbit, where a parabolic orbit has e = 1. 
/ Pericentric distance of the initial orbit 

pair time-scales, but is unlikely to affect the morphological 
disturbance time-scales which peak at the first pass and final 
merger (see §4.6 for discussion). 

The simulations only include feedback from supernovae. 
Most of these Sbc simulations were run with the 'stiff' n = 2 
supernova feedback equation of state; the Sbc prograde- 
prograde merger and radial orbit merger were run with the 
isothermal n = supernova feedback equation of state as 
well (SbcPPn=0, SbcRn=0). To test the simulations for con- 
vergence, the Sbc prograde-prograde merger was also run 
with four and ten times as many particles as our typical sim- 
ulations (SbcPPx4, SbcPPxlO). All of the G-G simulations 
were run with both supernova feedback models. Although 
feedback from an active galactic nucleus (AGN) may be im- 
portant for the properties of the merger remnants, such feed- 
back is not expected to affect the large-scale merger proper- 
ties and morphology until after the coalescence of the nuclei. 
Such feedback is driven by a rapidly accreting AGN which 
may appear as an optically luminous quasar during the end 
stages of the merger. Incorporating a quasar into the ra- 
diative transfer code poses a number of technical challenges 
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and is beyond the scope of this paper, but will be included 
in a future version of SUNRISE. As we discuss in §5, the 
exclusion of AGN feedback effects will not aifect the mor- 
phological disturbance time-scales calculated here. 



3 ANALYSIS 

We replicate the observations and measurements of real 
galaxy mergers as closely as possible. Our simulations pro- 
duce multi-wavelength images of galaxy mergers including 
the effects of dust and star-formation. Most current morpho- 
logical measurements of the evolution of the galaxy merger 
rate are done in the rest-frame B or SDSS g (e.g. Conselice 
ct al. 2003, Lotz et al. 2008). This is because high spatial 
resolution rest-frame 4000 — 5000 A imaging can be done 
for local galaxies from the ground, for 2 ~ 1 galaxies with 
the Hubble Space Telescope Advanced Camera for Surveys, 
and for z < 3 galaxies with the HST Near-Infrared Cam- 
era and Wide Field Camera 3. Therefore wc focus on SDSS 
g morphologies for purposes of this paper, as these simula- 
tions can be used to calibrate the morphologies of galaxies 
currently observed at < 2: < 3. In the following section we 
describe how the simulated SDSS g images are degraded and 
analysed to match real galaxy morphology measurements. 

3.1 Image degradation 

The SDSS g images are produced by SUNRISE for each sim- 
ulation for 11 isotropically positioned viewpoints as a func- 
tion of time from ~ 0.5 Gyr prior to the first pass to > 1 Gyr 
after the final coalescence in ~ 30 — 250 Myr timesteps de- 
pending on the merger state. The field of view of the output 
images ranges from 200 kpc during the initial stages and pe- 
riod of maximal separation to 100 kpc during the first pass, 
final merger and post-merger stages. The intrinsic resolution 
of the output SUNRISE g-band images is 333 pc. 

The images output by SUNRISE have no background 
sky noise and no seeing effects, although they do have parti- 
cle noise and Monte Carlo Poisson noise. We degrade these 
images to simulate real data, but do not attempt to mimic a 
particular set of galgixy survey observations. The measured 
morphology is dependent on the observed spatial resolution 
and signal-to-noise to the extent that low spatial resolution 
(> 500 pc per resolution element) and low signal-to-noise 
(average S/N per pixel < 3-5) introduce biases in the mor- 
phology values. However, above these limits, measured mor- 
phologies are not dependent on spatial resolution or signal- 
to- noise (LPM04) . We re-bin the images to 105 pc per pixel 
and convolve the images with a Gaussian function with a 
FWHM = 400 pc. This was done to simulate the effect of 
seeing but maintain as high spatial resolution as possible. 
The values where chosen to match the SDSS with 1.5" see- 
ing, 0.396" per pixel plate scale for a galaxy at a distance 
such that 1.5" ~ 400 pc. We also add random Poisson noise 
background to simulate sky noise but scale this noise to 
maintain a high S/N for the galaxies (> 20 per pixel within 
the Petrosian radius). Our choice to simulate images with 
spatial resolution < 400 pc per resolution element and S/N 
per pixel > 20 means that our results here can be generally 
compared to any rest-frame ~ 4700 A morphological study 
with sufficient spatial resolution and signal-to-noise. This 



corresponds to galaxies with a distance modulus < 35.0 ob- 
served from the ground with a seeing PSF FWHM ~ 0.8"or 
galaxies aX z < 0.25 observed with HST ACS and a PSF 
FWHM 0.14". In §4.2 we discuss how the results presented 
here apply to high-redshift galaxies observed with HST at 
spatial resolutions ^ 1 kpc, where the morphological biases 
from spatial resolution can be important. 



3.2 Morphology Measurements 

The degraded image for each snapshot and camera is treated 
as an independent observation with no prior information ex- 
cept for the central position of the two galaxies, which is used 
to track the galaxies' identities. Each image is run through 
an automated galaxy detection algorithm integrated into our 
IDL code. This algorithm is similar to but simpler than the 
detection and de-blending algorithm of SExtractor (Bertins 
& Arnouts 1996), which is optimized for large images with 
many objects but is ill-suited for images of one or two ob- 
jects. The image is smoothed by a ~ 4 kpc boxcar and ini- 
tial segmentation maps of galaxies are determined based on 
a fixed surface-brightness threshold equal to 2anoise- The 
number of objects detected in the 2a threshold map is com- 
pared to the number of objects detected in a 15ct thresh- 
old map. If more than 1 object larger than the smoothing 
length is detected in either map, a de-blending algorithm is 
applied. The largest objects detected in the high threshold 
map are grown using an image dilation algorithm (the IDL 
DILATE function) and a 5 x 5 pixel square shape operator, 
with a limiting surface brightness set to the 2cr threshold. 
The de-blending algorithm adopted here results in similar 
segmentation maps to those used in LPM04 for the sample 
of local mergers. 

If the centres of the merging galaxies are less than 10 
kpc apart, they are generally detected as a single object. 
If 2 distinct galaxies are detected, the detection segmenta- 
tion maps are used to mask out the other galaxy while each 
galaxy's morphology is measured. The output segmentation 
maps are visually inspected. Occasionally, the detection al- 
gorithm will assign a tidal dwarf galaxy as a second primary 
galaxy or fail to adequately mask out the secondary galaxy. 
In these cases, the masking is done by hand and the morphol- 
ogy code is re-run. For this paper, we ignore the properties 
of any tidal dwarfs produced in the merger and focus only 
on the merging galaxies and their remnants. 

The centres of each galaxy are estimated by minimizing 
the second-order moment of the pixels assigned to the detec- 
tion segmentation map, and ellipticities and position angles 
are determined using the IDL task FIT_ELLIPSE (Fanning 
2002). The projected separation Rproj is measured when two 
galaxies are detected. The initial guesses at the centre, el- 
lipticity, and position angle are then used to calculate the 
Petrosian radii in circular and elliptical apertures, concen- 
tration, asymmetry, dumpiness, the Gini coefficient, and the 
second-order moment of the brightest 20% of the light (see 
Lotz ct al. 2004 and Conselice 2003 for further details). 

The Petrosian radius is defined as the radius at which 
the ratio of the surface brightness at rp to the mean surface 
brightness within rp is equal to a fixed value, i.e. 
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where 77 is typically set to 0.2 (Petrosian 1976). Because the 
Petrosian radius is based on a curve of growth, it is largely 
insensitive to variations in the limiting surface brightness 
and S/N of the observations. 

Concentration is defined in slightly different ways by 
different authors, but the basic function measures the ratio 
of light within a circular inner aperture to the light within 
an outer aperture. We adopt the Bershady et al. (2000) def- 
inition as the ratio of the circular radii containing 20% and 
80% of the "total flux" : 



C = 5 1ogl0f^) 
\r2o/ 



(2) 



where rso and r2o are the circular apertures containing 80% 
and 20% of the total flux, respectively. For comparison to 
the most recent studies of galaxy concentration, we use Con- 

sclicc's (2003) definition of the total fiux as the fiux con- 
tained within 1.5 Vp of the galaxy's centre {as opposed to 
Bershady's definition as the flux contained within 2 rp). 
For the concentration measurement, the galaxy's centre is 
that determined by the asymmetry minimization (see be- 
low). Bulge-dominated early- type galaxies generally have 
high concentrations (C ~ 4-6), while late- type discs have 
low concentrations (C ~ 2-3). On-going mergers may show 
very low concentrations or very high concentrations depend- 
ing on the merger stage and the brightness of the central 
starburst. 

The asymmetry parameter A quantifies the degree to 
which the light of a galaxy is rotationally symmetric. A is 
measured by subtracting the galaxy image rotated by 180 
degrees from the original image (Abraham et al. 1994, Con- 
selice et al. 2000). 

- h8o{i,j)\ 



(3) 



where / is the galaxy's image and Jiso is the image rotated 
by 180 about the galaxy's central pixel, and -B180 is the 
average asymmetry of the background. A is summed over 
all pixels within 1.5 rp of the galaxy's centre. The central 
pixel is determined by minimizing A. The asymmetry due 
to the noise must be corrected for, and it is impossible to 
reliably measure the asymmetry for very low signal-to-noise 
ratio images (LPM04). Objects with very smooth elliptical 
light profiles have low asymmetries (A < 0.05). Galaxies 
with spiral arms are more asymmetric (A ~ 0.1-0.2), while 
extremely irregular and merging galaxies are often highly 
asymmetric (A > 0.35). 

The smoothness parameter S has been developed by 
Conselice (2003), inspired by the work of Takamiya (1999), 
in order to quantify the degree of small-scale structure. The 
galaxy image is smoothed by a boxcar of given width and 
then subtracted from the original image. The residual is a 
measure of the dumpiness due to features such as compact 
star clusters. In practice, the smoothing scalelength is chosen 
to be a fraction of the Petrosian radius. 



Bs 



(4) 



where Is is the galaxy's image smoothed by a boxcar of 
width 0.25 rp, and Bs is the average smoothness of the back- 
ground. Like A, S is summed over the pixels within 1.5 rp 
of the galaxy's centre. However, because the central regions 



of most galaxies are highly concentrated, the pixels within 
a circular aperture equal to the smoothing length 0.25 rp 
are excluded from the sum. S is correlated with recent star- 
formation (Conselice 2003). 

The Gini coefficient G is a statistic based on the Loronz 
curve, the rank-ordered cumulative distribution function of 
a population's wealth or, in this case, a galaxy's pixel values 
(Abraham et al. 2003). The Lorenz curve is defined as 



L(p) 



1 r 



F-\u)du 



(5) 



where p is the percentage of the poorest citizens or faintest 

pixels, F(x) is the cumulative distribution function, and X is 
the mean over all (pixel fiux) values Xi (Lorenz 1905). The 
Gini coefficient is the ratio of the area between the Lorenz 
curve and the curve of "uniform equality" (where L(p) = p) 
to the area under the curve of uniform equality (= 1/2). For 
a discrete population, the Gini coefficient is defined as the 
mean of the absolute difference between all Xf. 



G 



2Xn(: 



n n 

n — 1) 

EE 



Xi — JCi 



(6) 



i=i i=i 



where n is the number of people in a population or pixels 
in a galaxy. In a completely egalitarian society, G is zero, 
and if one individual has all the wealth, G is unity. A more 
efficient way to compute G is to first sort Xi into increasing 
order and calculate 



n 

-^-—^{2i-n-l)\Xi 
nin — 1 ^ — ' 



(7) 



(Glasser 1962). G is high for objects with very bright nuclei 
(G ~ 0.6), whether those galaxies are highly concentrated 
ellipticals or mergers with multiple bright nuclei. It is low 
for objects with more uniform surface brightnesses, such as 
late- type discs (G ~ 0.4). 

Because G is very sensitive to the ratio of low sur- 
face brightness to high surface brightness pixels, G must 
be measured within a well-defined segmentation map. For 
the purposes of measuring G and M20, we have chosen to 
assign pixels to the galaxy based on the surface-brightness 
at the Petrosian radius as measured in elliptical apertures 
(LPM04; see Abraham et al. 2007 for a similar approach). 
Note that the resulting segmentation map is not ellipti- 
cal, but rather traces the isophote that matches the mean 
surface-brightness at the elliptical Petrosian radius. The 
Petrosian radius is a reproducible quantity that is relatively 
independent of signal-to-noise and surface-brightness dim- 
ming effects (LPM04). Designating galaxy pixels based on 
signal-to-noise cuts (e.g. Law et al. 2007) will result in unreli- 
able G values. This is because galaxies with the same 'shape' 
or 'morphology' but difi'ering luminosities will have different 
measured G values when signal-to-noise is used to define 
the galaxy pixels. For example, a low surface-brightness ex- 
ponential disc will have fewer low surface-brightness pixels 
assigned to its segmentation map than a morphologically 
identical high surface-brightness exponential disc, resulting 
in a lower G for the low surface-brightness disc. Moreover, G 
values measured within segmentation maps based on signal- 
to-noise cuts are not repeatable because the measured G will 
depend as much on the noise properties of the observations 
as the intrinsic galaxy properties. The results presented in 
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this work are robust to these effects because the pixel maps 
used to measure G are based on the surface brightness at 
the Petrosian radius. 

The total second-order moment Mtot is the flux in each 
pixel, fi , multiplied by the squared distance to the centre of 
the galaxy, summed over all the galaxy pixels assigned by 
the segmentation map: 



Table 3. Merger Stages 



Mt. 



(8) 



where Xc,yc is the galaxy's centre. The centre is computed 
by finding Xc,yc such that Mtot is minimized. The second- 
order moment of the brightest regions of the galaxy traces 
the spatial distribution of any bright nuclei, bars, spiral 
arms, and off-centre star-clusters. A/20 is defined as the nor- 
malized second order moment of the brightest 20% of the 
galaxy's fiux. To compute M20, we rank-order the galaxy 
pixels by fiux, sum Mi over the brightest pixels until the 
sum of the brightest pixels equals 20% of the total galaxy 
fiux, and then normalize by Mtot'- 



M20 = log 



Mtot 



while 



, /> < 0.2 ftot 



(9) 



Here ftot is the total fiux of the galaxy pixels identified 
by the segmentation map and fi are the fluxes for each 
pixel i, ordered such that /i is the brightest pixel, /2 is the 
second brightest pixels, and so on. The normalization by 
Mtot removes the dependence on total galaxy fiux or size. 
M20 always has a value < 0. M20 is anti-correlated with C 
for normal galaxies, with low values for early- type galaxies 
(M20 < —2) and intermediate values for late-type galaxies 
(M20 ~ —1.5). Mergers with multiple nuclei have high M20 
values ( > —1). 

3.3 Definition of Merger Stages 

From the true three-dimensional separations of the galaxy 
nuclei, we determine the timestep of the closest approach 
during the first pass (tfp), maximal separation after the first 
pass (tmax), and final merger of the nuclei where St < 1 kpc 
(imerg). Thcsc times are given in Table [3]3] for each simula- 
tion. Based on these events, we deflne 6 merger stages; pre- 
merger, first pass, maximal separation, final merger, post- 
merger, and merger remnant. The 'pre-merger' stage is from 
t=0 to 0.5 tfp. The 'first pass' stage encompasses the first 
pass and starts at 0.5 tfp and ends at 0.5 {tfp + tmax). 
The 'maximal separation' stage starts at 0.5 {tfp + tmax) 
and ends at 0.5 {tmax + tmerg). The 'merger' stage starts 
at 0.5 {tmax + tmerg) aud cuds at tmerg+ 0.5 Gyr. The 
'post-merger' stage is defined as between tmerg+ 0.5 Gyr 
and tme.rg+ 1.0 Gyr, while the 'remnant' stage is at times 
more than 1 Gyr after the merger event (tmerg -1- 1.0 Gyr). 
We show composite SDSS u — r — z images for each of these 
stages for a prograde-prograde Sbc merger simulation as 
viewed face-on (Figure 1) and edge-on (Figure 2). 

3.4 Merger Classification and Time-Scales 

In Figure [31 we show the empirical criteria for merger classi- 
fication via G — M20 , G — A and C — A morphologies for lo- 
cal samples of visually classified galaxies (LPM04). All three 
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The time since the start of the simulation is given for each event 
that defines a particular merger stage (see §3.3 for definitions). 

morphological merger criteria are based on the Borne et al. 
(2000) observations of local ultra-luminous infrared galaxies 
(ULIRGs), a subset of which have been shown to be gas- 
rich mergers with mass ratios between 1:3 and 1:1 (Dasyra 
et al. 2006). LPM04 found that ULIRGs visually classified as 
mergers could be distinguished from the sequence of normal 
Hubble type galaxies with 



G > -0.115 M20 + 0.384 



G > -0.4 A + 0.66 or ^ > 0.4 



(10) 



(11) 



Asymmetry alone is also often used to classify merger candi- 
dates. The calibration of local mergers by Conselice (2003) 
finds the following merger criterion: 



A > 0.35 



(12) 



These are the merger criteria used to calculate the mor- 
phological observability time-scales throughout most of this 
paper. 

Galaxies at higher redshift cannot be imaged at as high 
spatial resolution as local galaxies even when observed with 
HST. The measured morphologies of galaxies at 2; > 0.25 
imaged with HST will have non- negligible biases as a re- 
sult of this lower spatial resolution (LPM04). Therefore the 
merger criteria have been adjusted to account for these bi- 
ases in HST data by Conselice et al. (2003), Conselice, 
Blackburne, & Papovich (2005), and Lotz et al. (2008). For 
z < 1.2 galaxies observed with HST ACS Wide-Field Cam- 
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Figure 3. G — M20 , G — A, and C — A relations for local galaxies measured by LPM04. Empirically, normal galaxies lie below the dashed 
lines in the G — M20 and G — A plots and to the left of the dashed line in the C — A plot. Mergers are identified as galaxies which lie 
above and/or to the right of these divisions. While 90% of local ULIRGs show visual signs of merger activity, ULIRGs with double nuclei 
show the strongest signatures in quantitative morphology. 
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Figure 4. Time v. Rproj, star-formation rate per object. A, G, M20, and C for the high-resolution prograde-prograde Sbc merger with 
stiff supernova feedback, including the effects of dust (SbcPPxlO). Each merger stage is marked with a different colour (pre-merger: red, 
first pass: green, maximal separation: blue, final merger: orange, post-merger: magenta, remnant: cyan). Open diamonds are for the one 
merging galaxy and the merger remnant; crosses are for the other merging galaxy. The mergers show strong morphological disturbances 
and peaks in the star-formation rate at the first pass (green points) and final merger (orange points). 



era at rest-frame --^ 4000A, the revised G — M20 merger clas- 
sification is: 

G > -0.14 M20 + 0.33 (13) 

(Lotz et al. 2008). Conselice et al. (2005) find that decreased 
spatial resolution and surface brightness dimming at z > 0.5 
can lower the measured A in irregular galaxies by 0.05-0.15 



(also Shi 2008). In §4.2, we find a median offset of -0.05 for 
A when our simulations are convolved to match the spatial 
resolution of HST ACS WFC F814W^ observations at z ~ 
1. Therefore we suggest a revised merger criterion for G — 
A and A for HST observations of high-redshift galaxies as 
follows: 

G > -0.4 A 4- 0.68 or ^ > 0.35 (14) 
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and 

A > 0.30 



0.5 Nc p{merg) 



(19) 



(15) 



We will use these merger criteria in §4.2 for the simulations 
convolved to match HST ACS observations of galaxies at 
~ 1. 

Close kinematic pairs are also probable merging sys- 
tems. Recent studies of local kinematic pairs have selected 
objects with 5 < Rp,oj < 20h~^ kpc (Patton et al. 2002, de 
Propris et al. 2005) and relative velocities < 500 km s~^, 
while studies of pairs out to 2 ~ 1.4 has chosen objects with 
10 < Rproj < 30, 50, and 100 h^^ kpc and relative veloci- 
ties < 500 km s"^ (Lin et al. 2004). We assume ft = 0.7 and 
we estimate the time-scales during which merging galaxies 
can be found as separate objects within 5 < Rproj < 20, 10 
< Rproj < 30, 10 < Rproj < 50, and 10 < Rproj < 100 /i"^ 
kpc. The simulated merging galaxies always have relative 
velocities < 500 km s~^. 

The galaxy merger rate F is defined as the number of on- 
going merger events per unit volume (pmerg divided by the 
time Tmerg for the merger to occur from the initial encounter 
to the final coalescence: 



r = 



(16) 



However, the number density of galaxies identified morpho- 
logically as galaxy mergers <j)merg will depend on the time- 
scale Tmerg duriug which the merger can be observed given 
the method used to identify it such that 



T 



'T 



(17) 



Therefore the galaxy merger rate F can be calculated from 
the observed number density of galaxy merger candidates 



as follows; 



4^merg Tmerg 
Tmerg Tmerg 



Tmerc 



(18) 



The effective observability time-scale T given in Equa- 
tion[T8]is a weighted average of the time-scales over all likely 
merger mass ratios, gas fractions, and orbital parameters. 
We do not calculate this global observability time-scale here 
because this may require a cosmological model for the dis- 
tribution of galaxy merger properties if the observability 
time-scale for a given method varies strongly. Instead, we 
present the first steps towards calculating the global observ- 
ability time-scales by exploring the the dependence of the 
time-scales on a wide range of galaxy merger properties. 

We calculate each simulation's average observability 
time-scale for the G — M20, G — A, and A criteria given above 
by averaging the results of the 11 isotropic viewpoints. Here- 
after we drop the prime notation and refer to the observ- 
ability time-scale for each simulation as T. Because we wish 
to determine the number density of merger events rather 
than the number of galaxies undergoing a merger, galax- 
ies that have not yet merged but identified morphologically 
as merger candidates are weighted accordingly. The time 
that each pre-merger galaxy is morphologically disturbed is 
summed (not averaged) to the time that the post-merger 
system appears disturbed. No such weighting is done for the 
close pair time-scales as this factor is generally included in 
the merger rate calculation (e.g. Patton et al. 2000): 



where 4> is the number density of galaxies within the mag- 
nitude range of the observed pairs, Nc is the average num- 
ber of companions within the observed magnitude range per 
galaxy, p{merg) is the probability that a galaxy pair will 
merge, Tpair is the time-scale for which merging galaxies 
will meet the close pair criteria, and 0.5 is the weighting 
factor that accounts for the double counting of pairs. 



4 RESULTS 

The equal-mass galaxy merger simulations span a wide range 
of physical properties. These include the relative orienta- 
tions and orbital parameters of the merging galaxies, the gas 
fraction and scalelength, the assumed supernova feedback 
prescription, and masses. In this section, we explore how 
important these physical parameters, the inclusion of dust, 
and simulation resolution are to the predicted morphologies. 
We present the viewing-angle averaged observability time- 
scales for G — M20, G — A, and A morphologies (Table [G]! 
and close pair projected separations Rproj < 20, 30, 50, and 
100 kpc for each simulation (Table [7]). Finally, the mor- 
phologies of the remnants observed 1 Gyr after the merger 
are calculated (Table 

For each simulation, we examine the projected sepa- 
rations Rproj, measured morphologies (G, M20, C, A, and 
5*), and star-formation rate per object as a function of time 
and merger stage for each simulation for all 11 viewpoints. 
The initial segmentation maps computed to identify each 
galaxy are used to compute the total star-formation rate 
for each galaxy at each timestep and camera. In general we 
only examine the morphologies of the output images that in- 
clude the effects of dust absorption and scattering. In Figure 
131 we show the evolution with time for the highest resolu- 
tion prograde-prograde Sbc merger simulation, SbcPPxlO. 
The merger stages are colour-coded, with red for pre-merger, 
green for first pass, blue for maximal separation, orange for 
merger, magenta for post-merger, and cyan for the remnant. 
The scatter at each timestep reflects the scatter in morphol- 
ogy with viewing angle. Before the final merger, the mor- 
phologies are measured separately for each galaxy (open di- 
amonds and crosses) . After the merger, the system is treated 
as a single object (open diamonds). 

Most of the parabolic-orbit mergers show trends of mor- 
phology and star-formation rate with merger stage similar 
to the SbcPPxlO merger (Figures 2] and [S]). Leading up to 
and including the first pass (green points) , the morphologies 
become increasingly disturbed as tidal tails are formed and 
the galaxies overlap in projection. The star-formation rate 
per object peaks at the first pass because the system appears 
as one object, but remains enhanced above the initial star- 
formation rates as the galaxies approach their maximal sep- 
arations (blue points). Dust starts to obscure star- formation 
in the nuclei, lowering G and G and increasing M20. Strong 
morphological disturbances are observed again at the final 
merger (orange points). The star- formation rate per object 
reaches its peak at or just after the final merger and gen- 
erally continues at high levels until ~ 0.5-1 Gyr after the 
coalescence of the nuclei. In Figure [5] we show the progres- 
sion of the SbcPPxlO merger in G - M20, G ~ A, and G - A 
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Figure 5. Top: Morphological evolutionary tracks in G — M20, G — A, and C — A space for the high resolution prograde-prograde 
Sbc merger (SbcPPxlO) averaged over all viewing angles. The initial galaxy (red), first pass (green), maximal separation (blue), final 
merger (orange), and last computed merger remnant (cyan) morphologies are labeled as 1-5 respectively. The black dashed lines show 
the empirical merger classification criteria from Figure 3. Bottom: The morphologies for all the timesteps and cameras for the SbcPPxlO 
simulation. The merger stages are indicted by the different colours, as in Figure |4] For this simulation, G — M20 identifies SbcPPxlO 
mergers primarily at the first pass, while G — A and A find mergers at both the first pass and final merger. 



averaged over all viewing angles (top) and for all 11 cam- 
eras (bottom). The initial galaxies (red points) start with 
late-type disc morphologies in all three plots. The system 
become disturbed in G — M20 space during the first pass 
(green points) and disturbed in G — A and C — A during 
the first pass and final merger (orange points). The post- 
mergers and remnants gradually end up with early-type disc 
morphologies in all three plots (cyan points). 

4.1 Numerical Resolution 

Most of our simulations were run with ~ 10^ particles 
per galaxy, with gravitational softening lengths of 400 pc 
and 100 pc for the dark matter and baryonic particles re- 
spectively (Table 1). The number of simulation particles 
affects both the spatial resolution of the simulation and 
the noise in the fluctuations of the gravitational potential. 
This number of particles was found to adequately recover 
the system-averaged star-formation histories and remnant 
properties when compared to simulations with larger num- 
bers of particles (Cox et al. 2006). Here we compare the 
time-dependent morphologies of the stiff supernova feed- 
back prograde-prograde Sbc merger (SbcPP) to simulations 
run with 4 and 10 times as many particles (SbcPPx4 and 
SbcPPxlO) to determine if the standard numerical resolu- 



Table 4. Morphological Timcscales v. Resolution 



Simulation 


T(G - M20) 
(Gyr) 


T(G - A) 
(Gyr) 


T(A) 
(Gyr) 


No corrections 


SbcPPxlO 

SbcPPx4 

SbcPP 


0.27± 0.13 
0.30± 0.12 
0.77± 0.15 


0.95± 0.16 
0.79± 0.17 
2.12± 0.51 


l.Olit 0.14 
0.82± 0.17 
2.74± 0.30 




t > 0. 


6 Gyr 




SbcPPxlO 

SbcPPx4 

SbcPP 


0.26± 0.10 
0.30± 0.11 
0.56± 0.19 


0.90± 0.14 
0.78± 0.17 
1.45± 0.31 


0.94± 0.13 
O.SOit 0.15 
1.94± 0.28 


t > o.e 


Gyr; 5M20 = 


-0.157; 5A = - 


0.115 


SbcPP 


0.39± 0.16 


0.78± 0.21 


0.74± 0.17 


PP=prograde-prograde; xlO 


= 10 Npart; x4 


= 4 Npart 



tion of the simulations is adequate also for analysing mor- 
phology. These higher numerical resolution simulations were 
processed by SUNRISE and the morphologies of the output 
SDSS g images were compared to the fiducial simulation 
at each timestep. The time evolution is slightly different 
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Figure 6. A Morphology v. Time for the prograde-prograde Sbc simulation run with lOx (SbcPPxlO) and Ix (SbcPP) the standard 
number of particles. The red error-bars show the standard deviation of the morphology differences within 0.25 Gyr bins. Prior to 0.6 
Gyr (dashed line), the low-resolution simulations has higher M20 and A values than the high-resolution simulation. 



in distinct simulation runs, so when the morphologies are 
rapidly changing, the morphologies, half-light radii (J?i/2), 
and galaxy separations may be significantly different for the 
different runs at a given timestep. However, for the major- 
ity of the simulation timesteps, differences in the morpholo- 
gies will reflect the differing spatial resolution and the noise 
in the gravitational potentials of the simulations. We find 
no resolution dependence for the time-scales of close pair 
projected separations Rproj, and so we focus only on the 
morphologies in this section. 

In Figure |6l we plot the difference between the stan- 
dard resolution simulation SbcPP and the highest resolu- 
tion simulation SbcPPxlO as a function of time for all 11 
viewing angles including the effects of dust. We find strong 
offsets in the half-light radii {R1/2) and morphologies for a 
few timesteps immediately after the first pass at 0.6 Gyr 
and immediately before the final merger at 1.7 Gyr, as ex- 
pected from small timing differences between the different 
simulations. For the majority of timesteps, the mean differ- 
ences between G, C, and effective radii for the standard and 
high resolution simulations are negligible but show signifi- 
cant scatter with viewing angle after the first pass. A/20, A 
and S do show systematic offsets particularly for the initial 
undisturbed galaxies (t < 0.6 Gyr; dashed vertical line in 
Figure [6]) . 

We compare the location of the brightest 20% of the 
pixels for the standard resolution and high resolution images 
of the initial discs in Figure [T] Young star particles in the 
spiral arms of the initial galaxies dominate the morpholo- 



gies because they are not adequately sampled in the fiducial 
simulations. Because star formation is implemented by cre- 
ating star particles whose mass depends on the resolution 
of the simulation, the fiducial simulations have fewer and 
brighter star clusters. For isolated and pre-merger galaxies, 
most new star particle formation happens at star-formation 
surface densities close to resolution of the output images 
(EsFfl ~ 0.03 M© yr"^ kpc"^, for 10^ Mq new star parti- 
cles, images with a spatial resolution ~ 400 pc, and 0/B 
star lifetimes ~ 20 Myr; see also Figure 3 in Cox et al. 
2006). Therefore the new stars are concentrated into a sin- 
gle particle within a single spatial resolution element, and 
the stochastic formation of individual star particles has a 
strong infiuence on the morphologies. As the numerical res- 
olution is increased and the mass of a new star particle de- 
creases, new stars are distributed over several particles. This 
decreases the typical surface brightness of the star-forming 
regions and the morphologies are less dominated by stochas- 
tic star particle formation. During the merger and remnant 
stages, most of the star-formation happens at gas densities 
well above the numerical limit of the fiducial simulations 
(see Figure 8 in Cox et al. 2006). Therefore the new stars 
are already distributed over multiple particles for the fidu- 
cial simulations, and the morphologies are not dominated 
by stochastic new star particle formation (Figure [Tj. 

After the first pass, the standard resolution simulation 
continues to show small but significant offsets to higher M20 , 
A, and S relative to the high resolution simulation (Fig- 
ure |6]). Even when we ignore the t < 0.6 Gyr timesteps, we 
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Figure 7. Top: SDSS g-band images for standard resolution sim- 
ulation SbcPP at time = 0.44 Gyr and 1.76 Gyr. Bottom: Same 
for high numerical resolution SbcPPxlO simulation. Thick eon- 
tours show the segmentation maps used to compute G and M20 ■ 
The thin contours show the pixels containing brightest 20% flux. 
Prior to the merger, the relatively low star-formation surface den- 
sities result in new stars that are concentrated in a single particle 
for a given star-forming region in the standard resolution simula- 
tions, resulting in artificially high M20, A, and S values. During 
and after the merger, high star-formation surface densities pro- 
duce multiple new star particles per star-forming region, resulting 
in more consistent morphologies between the standard and high 
numerical resolution simulations. 



still find morphological disturbance time-scales twice as long 
for the standard resolution simulations (Table |4}. We com- 
pute the mean offsets in M20 (— 0.157)and A (—0.115) be- 
tween the standard and high-resolution simulations after the 
first 0.6 Gyr. We recompute the standard-resolution merger 
time-scales ignoring the pre-merger initial galaxy morpholo- 
gies and correcting the morphologies by the mean offsets 
at t > 0.6 Gyr. These corrected time-scales are within la 
of the high resolution simulation merger time-scales, where 
Icr is the standard deviation derived from the 11 viewing 
angles (Table |4]). Given the large S offsets with simulation 
resolution, we do not include S in our analysis. 

We compare the output morphologies of the 4x and lOx 
resolution simulations (SbcPPx4, SbcPPxlO) to check that 
these simulations are resolved. Here the morphologies and 
time-scales agree quite well, even during pre-merger stage 
(Figure [HI Table [4]). This work is a first attempt at com- 
puting the wavelength-dependent morphologies for a large 
parameter space upon which future studies can build. There- 
fore we choose to examine only timesteps after the pre- 
merger stage [t > 0.6 Gyr) and apply the same corrections 
for M20 and A when calculating the observability time-scales 
for the standard-resolution simulations throughout this pa- 
per. 



4.2 Image Resolution 

At z > 0.25, even galaxies observed with HST have im- 
ages with worse spatial resolution than our fiducial resolu- 
tion (400 pc per resolution element). As we discuss in §3 
and show in LPM04, morphologies measured in images with 
spatial resolutions worse than 500 pc per resolution element 
have resolution-dependent biases. In principle, for images 
with low spatial resolution one should model the redshift and 
PSF-dependent biases for one's particular dataset. However, 
the turnover in the angular-size/ redshift relation is such 
that the angular scale of galaxies does not change dramati- 
cally at z > 0.6. In order to determine if the time-scales cal- 
culated here can be applied to HST images of high-redshift 
galaxies, we measure the morphologies and time-scales of 
the high numerical resolution simulation SbcPPxlO tuned to 
match typical HST observations of galaxies at z ~ 1. We use 
the TinyTinfl software to calculate the PSF of ACS WFC 
in the F9>14W (wide I) filter. We scale the PSF FWHM 
(0.14") to 1.1 kpc to match the angular scale at 2: ~ 1 and 
convolve the SbcPPxlO simulation images with this PSF. 

The morphologies measured from these images show 
non-negligible offsets from the images convolved to 400 pc 
resolution for M20, C, and A. The median offsets calculated 
for all timesteps and viewing angles are 5G = —0.01 ± 0.06, 
SAha = -0.03 ± 0.2, 5C = -HO. 12 ± 0.34 and 5A = -0.05 ± 
0.06. These offsets are consistent with the artificial redshift 
tests of real galaxies by LPM04 and Conselice et al. (2005). 
Without any corrections to the merger criteria given in Eqns. 
10-12, these shifts result in shorter observability time-scales 
(T(G - M20) = 0.14 ± 0.12 Gyr, T(G - A) = 0.58 ± 0.17 
Gyr, T(A) = 0.47 ± 0.17 Gyr vs. 0.26 ± 0.10, 0.90 ± 0.14, 
and 0.94 ± 0.13 Gyr, respectively). 

If revised merger criteria of Eqns. 13-15 (which take into 
account of the effect of decreased spatial resolution) are ap- 
plied, then the derived time-scales are in better agreement 
with the higher spatial resolution images: T(G — M20) = 
0.25 ± 0.12 Gyr, T{G - A) = 0.70 ± 0.17 Gyr, and T{A) = 
0.58 ± 0.15 Gyr. The pair time-scales are also ~ 200 Myr 
shorter than the fiducial resolution images. Although merg- 
ing objects may be more difficult to detect at high redshift 
with HST observations, the ~ 200-300 Myr offsets in the 
asymmetry and pair time-scales are significantly less than 
the variations associated with the gas properties and merger 
orbital properties (§4.3 and 4.4). We conclude that the time- 
scales presented here can be applied to HST observations 
of high-redshift galaxies without introducing uncertainties 
larger than those from the unknown distribution of merger 
properties. 



4.3 Dust and Viewing Angle 

The presence and distribution of dust has a strong effect on 
the measured morphologies starting at the first pass until 
after the final merger. In Figure [9] we plot the difference 
in the morphologies when dust is and is not included for 
SbcPPxlO. During the merger, dust mitigates the effect of 
star-formation on the morphologies. The presence of dust 
lowers G and C and increases M20 because the brightest 



^ J. Krist & R. Hook; ^ttp://www.stsci.edu/hst/software/tin3d:im| 
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Figure 8. A Morphology v. Time for the prograde-prograde Sbc simulation run with lOx (SbcPPxlO) and 4x (SbcPPx4) the standard 
number of particles. The red error-bars show the standard deviation of the morphology differences within 0.25 Gyr bins. The morphologies 
of the initial galaxies agree well at all times, including the initial stages at t < 0.6 Gyr (dashed line). 



Table 5. 


High Rcsolut 


on Prograde-Prograde Sbc Merger (SbcPPxlO) Timcscales v. Viewing Angle 




Cam 


T{G - A/20) 


T(G - A) 


T(A) 


T(5 < Rp < 20) 


T(10 <Rp< 30) 


T(10 <Rp< 50) 


T(10 <Rp< 100) 




(Gyr) 


(Gyr) 


(Gyr) 


(Gyr) 


(Gyr) 


(Gyr) 


(Gyr) 





0.24 


0.88 


0.98 


0.05 


0.24 


0.93 


1.30 


1 


0.32 


0.86 


0.86 


0.00 


0.22 


1.05 


1.20 


2 


0.44 


1.15 


1.15 


0.49 


0.29 


0.46 


0.61 


3 


0.27 


0.76 


0.98 


0.00 


0.10 


0.68 


1.20 


4 


0.10 


0.68 


0.71 


0.27 


0.73 


1.03 


1.03 


5 


0.39 


0.95 


0.86 


0.17 


0.32 


0.88 


1.20 


6 


0.24 


1.15 


1.05 


0.10 


0.24 


0.93 


1.34 


7 


0.27 


0.83 


1.10 


0.15 


0.83 


0.93 


1.25 


8 


0.15 


0.86 


0.81 


0.05 


0.29 


1.08 


1.25 


9 


0.20 


0.86 


0.86 


0.10 


0.36 


1.10 


1.25 


10 


0.20 


0.90 


1.00 


0.00 


0.22 


0.88 


1.05 


mean 


0.26± 0.10 


0.90± 0.14 


0.94± 0.13 


0.15 ± 0.18 


0.39 ± 0.22 


1.00 ± 0.17 


1.22 ± 0.20 



Computed for timesteps t > 0.6 Gyr, as discussed in §4.1. The camera angles are given in §2.2. 



star-forming regions near the centres of the merging galaxies 
are enshrouded. The G—M20 observabiUty time-scale is most 
strongly affected by the presence of dust, and is a factor of 2 
less when dust is included ( Table The G — A and A time- 
scales are relatively independent of extinction because the 
measured asymmetry is less affected (Figure O . The close 
pair time-scales are unchanged, as the measured positions 
and projected separations are unaffected by dust extinction. 
Most of the gas-rich simulations presented here continue to 



form stars at a significant rate (> 2 Mq yr"'^) at 1 Gyr after 
the merger of the nuclei, and retain significant amounts of 
gas and dust. If the dust is ignored, the remnants appear 
highly concentrated and relatively blue (TablelS)). With dust, 
the remnants appear less concentrated with lower C, G, and 
higher M20 and A values because the central star-formation 
is obscured and dust increases the asymmetry of the merger 
remnant (Table [8)|. 

The morphologies observed in SDSS g-band depend on 
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Figure 9. A Morphology v. Time for the high numerical resolution progrado-progrado Sbc merger (SbcPPxlO) with and without the 
effects of dust extinction. 
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Figure 10. Morphology as a function of viewing angle for the 
high resolution prograde-prograde Sbc merger (SbcPPxlO) . The 
red symbols show the initial galaxy morphologies and the cyan 
symbols show the final remnant morphologies. The measured 
morphologies of the initial disc galaxies do not change signifi- 
cantly with viewing angle. The merger remnant shows signifi- 
cantly higher M20 and A values when the dust lane is viewed 
edge-on (cameras 3 and 4). 



the viewing angle, in part because the dust lanes will pref- 
erentially absorb blue hght along certain lines of sight. Pro- 
jection effects and the relative orientations of the merging 
galaxies will also change the projected separations and ob- 
served morphology. The scatter in the morphology at a given 
timestep in Figure|4]is the result of the different viewing an- 
gles for the 11 different SUNRISE cameras. This scatter is 
largest immediately after the first pass and during the final 
merger when the system is most asymmetric, and is smallest 
for the merger remnant which is more spherically symmet- 
ric. We show the dependence of the measured morphologies 
on viewing angle in Figure [10] for the initial disc galaxies 
and the final remnant. The measured morphologies do not 
change significantly with viewing angle for the initial galax- 
ies. The remnant shows significantly higher M20, A, and 
S values when the final dust lane and star-forming disc is 
viewed edge-on (cameras 3 and 4) . In Table (5] we give the 
SbcPPxlO simulation merger time-scales for each viewing 
angle including the effect of dust. The standard deviation 
over all 11 viewing angles is ~ 100 Myr for T{G — M20), 
T{G — A), and T{A). The close pair time-scales and pro- 
jected separations Rproj also depend on viewing angle, as 
the galaxies will have smaller Rproj along some lines of sight. 
The close pair time-scales have a standard deviation ~ 200 
Myr over the 11 different viewing angles. 



4.4 Orientation and orbital parameters 

We examine the importance of the orbits and relative ori- 
entations of the merging galaxies to their morphologies. 
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Figure 11. Time v. Rprojj star-formation rate per object, A, G, M20, and C for the standard resolution prograde-prograde Sbc merger 
(SbcPP), the prograde-retrograde Sbc merger (SbcPR), the retrograde-retrograde Sbc merger (SbcRR), the polar orientation Sbc merger 
(SbcPol), and the radial orbit Sbc merger (SbcR). Each merger stage is colour-coded as in Figure |4] The parabolic orbit Sbc mergers 
show similar peaks in star-formation and morphological disturbances at the first pass and final merger, while the radial orbit Sbc merger 
(SbcR) has a single peak in both star-formation and disturbed morphology after the initial encounter. 



Sbc mergers initialized on parabolic orbits with pericen- 
tric distances Rperi = 11 kpc and stiff feedback were simu- 
lated with roughly prograde-prograde (SbcPP), prograde- 
retrograde (SbcPR), retrograde-retrograde (SbcRR), and 
polar (SbcPol) orientations. All of these simulations have 
similar orbital decay times, with the polar orientation 
merger taking a few 100 Myr longer for the final merger 
to occur (Table 3). All of the parabolic Sbc simulations 
experience peaks in star-formation and asymmetry at the 
first pass and the final merger, with the maximum star- 
formation rate depending on the relative orientation of the 



discs (Figure [TT|) . The strength of the morphological distur- 
bances also depends on the orientation of the galaxies, with 
the intrinsically asymmetric polar and retrograde-prograde 
mergers showing the highest asymmetries (Fig. Ilip . The 
viewing-angle averaged time-scales during which a partic- 
ular set of quantitative morphologies are disturbed vary by 
a factor of 2 for these different orientations {T{G ~ AI20) 
~ 0.3-0.6 Gyr, T{G - A) ~ 0.8-1.3 Gyr, and T{A) ~ 0.7- 
1.5 Gyr; Table [S] Figure [T^ . This is in reasonable agree- 
ment with the ~ 0.7-1.0 Gyr asymmetry time-scales of the 
star-particles of equal-mass prograde-inclined, retrograde- 
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Figure 12. G - A/20, G ~ A, and C - A for the same simulations as Figure [TTI fSbcPP. SbcPR, SbcRR, SbcPoI, and SbcR). Each 
merger stage is colour-coded as in Figure |4] The SbcRR and SbcR simulations are more likely to have disturbed morphologies during 
the maximal separation stage between the first pass and final merger (blue points) than the other Sbc simulations. 
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Figure 13. Simulated merger remnant morphologies, including the effects of dust, for G— M20 , G — A, and C ~ A. The Sbc-Sbc remnants 
are black asterisks, and the G-G remnants are black diamonds. Simulations run with n = supernova feedback are surrounded by a black 
square. Also plotted are the SDSS g or B morphologies of local galaxies from SDSS and the Frei et al. catalog, measured by LPM04. 
The remnants of all the simulations run with stiff n = 2 supernova feedback have quantitative morphologies similar to Sb galaxies (green 
triangles). Some simulations with isothermal n = supernova feedback approach E/SO morphologies (red diamonds). 



inclined, and prograde-retrograde merger simulations found 
by Conselice (2006). The retrograde-retrograde merger is 
disturbed for the longest time for all of the quantitative 
morphology measures. The time-scales are also sensitive to 
the criteria used to identify the merger (Table ID Figure fT2|) . 



The typical T{G - M20) is ~ 0.4 Gyr, while T{A) ~ 1.1 
Gyr and T{G — ^) ~ 1.0 Gyr. The close pair time-scales 
do not vary strongly with orientation {ST ~ 200 Myr; Ta- 
ble [7]) but do depend on the range of Rproj chosen, with 
Rproj < 30 kpc time-scales often significantly shorter than 
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Table 6. Equal-Mass Merger Morphological Timescales 



Simulation 


T(G-M2o) 


T{G-A) 


T{A) 




(Gyr) 


(Gyr) 


(Gyr) 




Sbc-Sbc mer^ 


5ers 




oDCJr^Jr^XiU \n\J QUSt ^ 


U.44it yj.io 


n QO_i_ n QA 


1 1 o_i_ n /I n 

l.lZlt U.4U 


SbcPPxlO 


0.26± 0.10 


0.90± 0.14 


0.94± 0.13 


SbcPPx4 


0.30± 0.11 


0.78± 0.17 


0.80± 0.15 


SbcPP 


0.39± 0.16 


0.78± 0.21 


0.74± 0.17 


SbcPR 


0.31± 0.10 


0.98± 0.11 


1.12± 0.11 


SbcRR 


0.60± 0.18 


1.33± 0.32 


1.46± 0.31 


SbcPol 


0.40± 0.25 


l.lOit 0.24 


l.lOit 0.29 


SbcPPr- 


0.57± 0.20 


0.77± 0.19 


0.70± 0.16 


SbcPPr-l- 


1.03± 0.74 


0.93± 0.48 


1.19± 0.57 


SbcR 


0.44± 0.21 


0.93± 0.30 


0.61± 0.12 


SbcPPn=0 


0.42± 0.20 


0.90± 0.35 


0.76± 0.17 


SbcRn=0 


0.23± 0.09 


0.86± 0.35 


0.42± 0.09 


G-G mergers 


G3PP 


0.17± 0.06 


0.32± 0.07 


0.22± 0.11 


G2PP 


0.22± 0.14 


0.31± 0.15 


0.25± 0.19 


GIPP 


0.24± 0.05 


0.35± 0.14 


0.30± 0.16 


GOPP 


0.30± 0.09 


0.43± 0.17 


0.39± 0.16 


G3PPn=0 


0.19± 0.08 


0.38± 0.15 


0.19± 0.16 


G2PPn=0 


0.22± 0.06 


0.30± 0.12 


0.21± 0.10 


GlPPn=0 


0.61± 0.24 


0.45± 0.13 


0.30± 0.12 


GOPPn=0 


0.31± 0.09 


0.48± 0.14 


0.41± 0.17 



Computed for timesteps t > 0.6 Gyr, as discussed in §4.1. The 
standard resolution simulations also have SM20 = —0.157, 
5A = —0.115 correction applied, as in the last row of Table|4] 

typical T{G — A) and T{A) values. All of the simulations 
show enhanced star-formation for significantly longer peri- 
ods than the morphological disturbances, with the peaks 
in the star-formation rates often occurring after the peaks 
in the asymmetry. The remnant morphologies are generally 
similar (Table|8}, although the retrograde-retrograde merger 
remnant has a higher A/20 value than the prograde-prograde 
merger ( — 1.77 v. —1.93) and both the retrograde-retrograde 
and polar merger remnants have higher asymmetries (0.15 
V. 0.0). The remnant morphologies are more consistent with 
early-type disc galaxies (Sb) than spheroids (E/SO) (Figure 
I13p . The remnants are forming stars at ~ 5 — 6 Mq yr~^. 

A sub-parabolic Sbc-Sbc merger simulation with a 
highly radial orbit, zero net angular momentum, prograde- 
retrograde orientation and stiff feedback was also examined 
(SbcR). Because the galaxies start with low relative veloc- 
ities, it takes significantly longer for the first pass to oc- 
cur (1.3 Gyr v. 0.6 Gyr). However, there is significantly 
less time between the first pass and the final merger (0.3 
Gyr V. 1.1 Gyr; Table[S}. As a result, the morphologies and 
star-formation rates as a function of merger stage are quite 
different from the parabolic orbit simulations (Fig. 1111 112|) . 
Asymmetry peaks during the first pass and maximal sepa- 
ration stages, but is quite low during the final merger. On 
the other hand, the star-formation rate is strongly enhanced 
throughout the later merger stages and reaches its peak dur- 
ing the final merger. The G — M20, G~ A, and A time-scales 
are similar to the parabolic -Rperi = 11 Sbc simulations. 
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Figure 14. Time v. Rproji star-formation rate per object. A, G, 
M20, and C for the standard resolution prograde-prograde Sbc 
mergers with small Rperi (SbcPPr-) and large Rperi (SbcPPr-|-). 
Each merger stage is colour-coded as in Figure |4] The large Rperi 
simulation takes significantly longer to merge, and has lower peak 
star-formation rates and asymmetries. 



The close pair time-scales, however, are naturally ~ 50% 
shorter than parabolic orbits. The merger remnant has a 
large bulge surrounding by a very blue star-forming ring. 
It has the highest star- formation rate (13 M© yr~^) of any 
of the simulations, and because of the bright blue ring, its 
morphology is the most asymmetric (0.25) and disc-like in 
its G, M20 and C values. 

Two additional prograde-prograde parabolic orbit simu- 
lations with smaller and larger pericentric distances {Rpsri ~ 
5.5, 44 kpc V. 11 kpc) were analysed. The small Rpei i sim- 
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ulation (SbcPPr-) takes 340 Myr less to merge, while the 
large Rperi (SbcPPr+) simulation takes an additional 2 Gyr 
to merge (Table |6l Figure I14|l . The merger time-scales and 
properties of the small Rperi simulation are similar to the 
fiducial SbcPP simulation, with somewhat shorter close pair 
time-scales for Rproj < 50 and 100 kpc (Table [T] Figure 
I14p . The large Rperi simulation, however, has little enhanced 
star-formation and lower asymmetries during the first pass, 
and experiences less enhanced star-formation during the fi- 
nal merger because more gas has been consumed in 'normal' 
disc star-formation (Fig. I14|l . Despite weaker morphological 
disturbances, the G — M20 observability time-scales for the 
large Rperi simulations is significantly larger (1.0 Gyr v. 0.4 
Gyr) but with larger scatter with viewing angle. The G — A 
and A time-scales are also longer with larger scatter, as are 
the the close pair time-scales for Rproj < 50 and 100 kpc. 
The large Rperi remnant has somewhat higher G, M20 and 
C values than the Rperi ~ 5.5 and 11 kpc remnants, mak- 
ing its morphology more like typical spheroidals (Table |8]). 
Visual inspection of the remnant shows that the recent star- 
formation in the large Rperi is much more centralized, while 
the smaller Rperi simulations have an extended disc of young 
stars. Because the SbcPPr-|- simulation takes twice as long 
to merge, most of its remnant's cold gas has been consumed 
in star-formation during the merger. 

In summary, we find that orientation and large pericen- 
tric distances can have a significant effect on the time-scales 
during which mergers can be identified morphologically. 
Some relative orientations of the merging system increase 
the strength of the morphological disturbances (prograde- 
retrograde, polar), while other orientations increase the 
time-scales of those disturbances (retrograde-retrograde). 
Large pericentric distances naturally result in long orbital 
decay times, which suggests that the duration of the merger 
is as important to the time-scales of morphological distur- 
bances as the orientation of the merging galaxies. The or- 
bits do affect the timing of the morphological disturbances. 
Most of the parabolic-orbit simulations show peaks in mor- 
phological disturbances at the first pass and final merger. 
The highly radial orbit simulation shows a single peak in 
asymmetry during the first pass, while the large pericen- 
tric radius simulation shows a less dramatic enhancement of 
asymmetry during the first pass and final merger. However, 
the highly radial orbit simulation also has morphological dis- 
turbance time-scales that agree with the parabolic orbits to 
within the scatter associated with the viewing angle. En- 
hanced star-formation rates generally occur for longer du- 
rations than the morphological disturbances, with the star- 
formation rates often peaking after the asymmetries. The 
remnants all have similar concentrations (G, M20, and C) 
consistent with early-type spirals but only the retrograde- 
retrograde, polar, and highly radial orbit merger remnants 
show significant asymmetries {A > 0.1). 

4.5 Gas fraction and scalelength 

Gas-rich mergers undergo significant star-bursts triggered 
by the varying tidal forces and infiow of gas during the 
merger. Because these new stars can influence the measured 
morphologies, it is likely that the amount of gas available to 
form stars affects morphology time-scales during the merger 
process. The Sbc galaxy and the G3 galaxy have similar total 
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Figure 15. Time v. Rproj, star-formation rate per object. A, 
G, M20, and C for the prograde-prograde 1.2 X lO^^ Mq G3 
merger (G3PP) and less massive prograde-prograde 5 X lO^'^ Mq 
G2 merger (G2PP). Each merger stage is colour-coded as in Fig- 
ure |4] The G3 merger is less morphologically disturbed than the 
Sbc prograde-prograde merger (SbcPP) during the first pass, and 
experiences a peak in the star-formation rate well after the peak 
in asymmetry at the final merger. 

masses (8.1x10^^ and 1.2x10^^^ respectively), and simi- 
lar bulge to disc stellar mass ratios (0.25 and 0.21). However, 
the Sbc galaxy has a much larger gas reservoir with over 50% 
of its baryons in gas. The G3 galaxy, on the other hand, has 
only ~ 20% of its baryons in gas. Both models assume that 
the gas disc has a scalelength three times the scalelength 
of the stellar disc, but the She's gas and stellar discs are 
twice the adopted scalelengths for the G model. The end 
result is that the Sbc merger simulations have much more 
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Figure 16. Time v. Rproji star-formation rate per object, A, G, 
M20, and C for the prograde-prograde low-mass 2 X 10'^'^ Mq G1 
and 5 X 10^° Mq GO mergers (GIPP, GOPP). Each merger stage 
is colour-coded as in Figure |4] The lower mass mergers undergo 
less star-formation but have time-scales for disturbed morphology 
similar to the more massive G3 and G2 mergers. 



gas at large radii as well as a higher density of gas within 
the central regions. Cox et al. (2008) found that higher cen- 
tral gas densities lead to less merger-induced star-formation 
when compared to the undisturbed disc star-formation. The 
SbcPP merger has less efficient merger-driven star-formation 
throughout the interaction relative to the G3PP merger 
(where efficiency refers to the fraction of total gas converted 
to stars), with the SbcPP and G3PP mergers showing 23% 
and 46% more star-formation than their undisturbed coun- 
terparts respectively. However, the SbcPP merger experi- 



ences higher star-formation rates in general and at the first 
pass and final merger in particular, refiecting its high gas 
fraction. 

We compare the time-dependent morphologies and pro- 
jected separations for prograde-prograde parabolic orbit 
mergers with stiff feedback for the Sbc and G3 galaxies 
(SbcPP and G3PP; Figs. [11] and [15|. We find that the 
merger time-scales and morphologies are also affected by 
the gas disc properties. Despite similar pericentric distances, 
the G3PP simulation takes about 700 Myr longer for the nu- 
clei to coalesce than the SbcPP simulation. As a result, the 
time-scales during which the merging galaxies can be iden- 
tified as a close pair are also longer. The G3PP simulation 
spends three times longer as a very close pair {Rproj < 30 
h^^ kpc) than the SbcPP simulation (Table [7|. However, 
the time-scales for morphological disturbances are shorter 
for the lower gas fraction G3PP simulation by a factor of ~ 
2-4 (Table[Bl). Although the G3PP merger is more efficient at 
turning the available gas into stars than the SbcPP merger, 
the G3PP merger has less star-formation overall because of 
its larger bulge and lower gas fraction (Fig. [TJl I17p . The 
G3PP merger experiences less star-formation along tidal 
arms and lower asymmetries at the first pass than the SbcPP 
merger. Nevertheless, the remnants have similar quantita- 
tive morphologies and star-formation rates (Table [8|. Both 
simulations produce dusty remnants with G, M20, C, and 
A values consistent with bulge-dominated spirals (Fig. I13p . 
Both remnants have significant residual star formation (3-5 
Mq yr~^) and dust reddening. 

4.6 Mass 

The total mass involved in the merger may also affect mor- 
phologies and star-formation rates, as larger galaxies have 
deeper potential wells and produce stronger tidal forces. 
Equal-mass prograde-prograde merger simulations spanning 
a factor of 23 in total mass and a factor of 50 in stellar 
mass were examined to explore the effects of merger mass 
(G3PP, G2PP, GIPP, and GOPP simulations). The progen- 
itor galaxies have increasing gas fractions and total mass to 
light ratios with decreasing mass (Table 1). Both supernova 
feedback models were explored for all the G-series merg- 
ers. The orbits for all of the G-series mergers are slightly 
sub-parabolic, with eccentricities e = 0.95. This significantly 
shortens the decay times for the Gl and GO mergers, and 
hence may result in shorter close pair time-scales than would 
be observed for e = 1.0 orbits. We do not expect this to im- 
pact the morphology observability time-scales, as disturbed 
morphologies are apparent only at the first pass and final 
merger stages. The initial separations are less than 100 
kpc for the G2, Gl, and GO mergers, and so close pair time- 
scales are not computed when the initial separation is less 
than the measured range of projected separations (Table [7]). 

We find that all of the equal-mass G-series mergers show 
similar correlations of the morphologies with merger stage. 
The morphological disturbance time-scales are ^ 100 -200 
Myr longer for the lowest mass merger (GOPP) than the 
highest mass merger (G3PP; Table[6l Fig. 1151 [16)1 . However, 
the time between the first pass and coalescence of the nuclei 
is ~ a factor of 2 less for the lower mass merger than the 
highest mass merger (Table [6] Figure [T7[) . and the close pair 
time-scales at 10 < Rproj < 30 kpc reflect this (Table 
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Figure 17. G — M20, G - A, and C - A for the prograde-prograde G3PP, G2PP, GIPP, and GOPP simulations. The simulations span 
a factor of 23 in total mass, where the virial mass of the initial galaxy is 1.2 X 10^^ Mq for G3, 5.1 X 10^^ Mq for G2, 2.0 X 10^^ A^q for 
Gl, and 5.1 x 10^^ for GO. Each merger stage is colour-coded as in Figure|4] Unlike the Sbc mergers, the G-series simulations are only 
detected at the final merger (orange points). 



[T]). The close pair time-scales at 5 < Rproj < 20 kpc 
show longer time-scales for the lowest mass merger. This is 
an artifact of the object detection algorithm. Larger galaxies 
are most likely to be counted as one object at small separa- 
tions because they have larger scalelengths and overlapping 
isophotes. The merger remnants for all but the lowest mass 
merger (GOPP) are remarkably similar in their morphology, 
with the lower mass remnants showing lower star-formation 
rates and extinctions. The GOPP remnants are more disc- 
like in their G and M20 values, even when n = feedback is 
adopted. 



4.7 Supernova feedback 

The parabolic prograde-prograde and radial prograde- 
retrograde Sbc simulations and all of the G-series simula- 



tions were run with both supernova feedback models. Al- 
though the total gas consumption and star formation during 
the merger are similar for both feedback scenarios (Cox et al. 
2006, 2008), the isothermal n = feedback simulations with 
parabolic orbits experience more star-formation during the 
first pass and have less gas available for a second starburst 
during the final merger (Cox et al. 2006). The opposite is 
true for the highly radial orbit Sbc simulations because there 
is not enough time between the first pass and final merger 
for the low feedback simulation to consume large amounts of 
gas. The n = radial-orbit Sbc merger experience the high- 
est peak star-formation rate of all the merger simulations 
during the final merger (> 500 M© yr~^), while the n — 2 
radial-orbit Sbc merger experiences its peak star-formation 
rate during the first pass (~ 80 M© yr~^). 

We find that close pair time-scales and morphological 
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Figure 18. A Morphology v. time for the simulations with different supernovae feedback prescriptions (SbcPP = prograde-prograde 
Sbc; SbcR = radial orbit Sbc; G3PP = prograde-prograde G3). The M20 values are higher and the concentrations are lower for n = 2 
'stiff' supernovae feedback simulations during and after the final merger stage (orange, magenta, and cyan points). Each merger stage is 
colour-coded as in Figure |4] 



disturbances as a function of merger stage and merger time- 
scales are generally similar for the different feedback models, 
given the scatter with viewing angle (Tables [S] [7J . The pri- 
mary difference between the n = 2 stiff feedback and the 
n = isothermal feedback simulations appears in the prop- 
erties of the remnants (TableO Fig. ll3|l . The n = feedback 
remnants are significantly more like E/SO in the quantitative 
morphologies (G, M20 and C) because they have ~ 40-50% 
lower gas metallicities and hence less dust to obscure the 
nuclei. The lower gas metallicities are probably an artifact 
of our chemical enrichment scheme, rather than a robust 
prediction of the dust evolution. Supernova are assumed to 



produce only metals which enrich surrounding gas particles, 
but do not produce any gas particles themselves. If all sur- 
rounding gas particles are consumed in star- formation, the 
metals which would have be produced in supernova have no 
place to go and remain locked up in the stars. The n = 
feedback models experience more intense star-formation and 
consume a great amount of their existing gas during the 
first pass, hence are more effected by these limitations in 
our model. Nevertheless, while it is unclear if the dust prop- 
erties of the merger remnants will be strongly affected by 
feedback, it is likely that the n = 2 feedback remnants have 
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too much dust, as we do not include any dust destruction 
mechanisms. 



5 DISCUSSION 

Every equal-mass gas-rich merger simulation presented here 
exhibits quantitatively disturbed morphologies at some 
point along the merger process. However it is clear that 
quantitative morphological classifications based on G, M20, 
and A are sensitive only during the first pass and final 
merger stages for gas-rich equal-mass mergers, and will miss 
many interacting galaxies observed between the first pass 
and final merger as well as many recently merged systems. 
This is in qualitative agreement with the G, M20, and A 
values and merger stages of local ULIRGs. Two-thirds of 
the z ~ 0.1 ULIRG sample used to calibrate the G — M20, 
G — A, and C — A diagrams in LPM04 exhibit double or 
multiple nuclei, and therefore are merging systems observed 
at final merger stage before the coalescence of their nuclei 
or immediately at the first pass when the galaxies appear 
overlapping in projection (Figure 3). The G — M20, G — A, 
and C — A merger classification cuts used in LPM04 and this 
work identify 93%, 80%, and 76% of the double and multiple 
nuclei ULIRGs respectively. The detection efficiency is sig- 
nificantly lower for the single nucleus ULIRGs (46%, 71%, 
and and 54% for G — M20, G - A, and C - A respectively). 
This is also in reasonable agreement with our results here, 
assuming that single nucleus ULIRGs are observed after the 
first pass or after the coalescence of the nuclei. 

The duration, strength, and timing of the observed mor- 
phological disturbances depend on the merger orientation 
and orbital parameters, the gas properties of the initial 
galaxies, and the presence of dust. When dust is included, 
the merger observability time-scales depend most strongly 
on the gas properties, pericentric distance, and relative ori- 
entation. Galaxies with high gas fractions have more star- 
formation along tidal arms, producing stronger asymmetries 
during the first pass. Mergers with large impact parameters 
have long orbital decay time-scales, and exhibit disturbed 
morphologies for longer. Retrograde-retrograde mergers also 
show disturbed morphologies for 50-100% longer than 
prograde-prograde and prograde-retrograde mergers. We 
find that the supernova feedback prescription and the to- 
tal mass of the merging galaxies do not have a strong ef- 
fect on the overall duration of morphological disturbances. 
The relative orientations affect the strength of the morpho- 
logical disturbances, with the prograde-retrograde and po- 
lar orientations showing the strongest disturbances. The or- 
bital parameters and gas fractions have the strongest infiu- 
ence on the timing of the morphological disturbances. Most 
of the high gas-fraction (Sbc) parabolic orbits show mor- 
phological disturbances at the first pass and final merger, 
while the high gas-fraction highly radial orbit and large peri- 
centric distance simulations have weak disturbances at the 
first pass and stronger disruptions at the final merger. The 
lower gas-fraction (G-series) parabolic orbit simulations ex- 
perience less star-formation and morphological disturbances 
during the first pass, and hence are most likely to be de- 
tected morphologically during the final merger. 

Obscuration from dust has a very strong impact on the 
measured morphologies throughout the merger process until 



at least 1 Gyr after the coalescence of the nuclei. Dust extinc- 
tion is highest for the central nuclei where the star-formation 
rates are highest. Because much of the central light is masked 
by dust, this results in lower G, C, and higher M20 values 
during and after the merger. Dust lanes in the remnants 
can also produce higher asymmetries. However, the inclu- 
sion of dust does not significantly change the morpholog- 
ical disturbance time-scales during the prograde-prograde 
Sbc merger. Our models may overestimate the dust content 
in the remnants as dust is not destroyed by shock-heating 
nor is gas removed in a post-merger 'blowout' by an active 
galactic nucleus as predicted by other galaxy merger mod- 
els (e.g. Hopkins et al. 2006). However, dust destruction is 
expected to be most important after the final merger, when 
the star-formation rate and dust production have sufficiently 
declined. Therefore, while the dust and gas content of our 
merger remnants may be overestimated, the morphologies 
and time-scales calculated during the merger are unlikely to 
be affected by the destruction/removal of dust at late stages. 

The observability time-scales clearly depend on the 
method used to select merger candidates. The time-scale 
during which a merging system is a close pair at a par- 
ticular projected separation is not the same as the time- 
scale during which the system shows high asymmetries or 
high G — M20 values. The gas-rich Sbc mergers have 2-4 
times longer T(G - A) and T{A) than T{G - M20) and 
T(10 < Rproj < SOh~^ kpc). The lower gas-fraction G-series 
simulations, ou the other hand, have similar time-scales for 
G — M20, G — A, and A disturbances. Although G and M20 
are the most affected by dust, the G — M20 time-scale is 
the least affected by the merger parameters with typical 
time-scales ^ 0.2 — 0.3 Gyr for the G-series simulations and 
^ 0.3 — 0.6 Gyr for the Sbc simulations. The typical G — A 
time-scales are ~ 0.3 — 0.4 Gyr for the G-series simulations 
and ~ 0.8 — 1.1 Gyr for the Sbc simulations. The typical 
A time-scales are ~ 0.2 — 0.3 Gyr for the G-series simula- 
tions and ~ 0.7 — 1.1 Gyr for the Sbc simulations. Thus the 
G — M20 time-scales have ~ 0.4 Gyr dispersion, while the 
G — A and A time-scales have ~ 0.8 Gyr dispersion. While 
more mergers may be identified using asymmetry given the 
longer asymmetry time-scales, the merger rate calculated us- 
ing G — M20 mergers will be less uncertain given the better 
consistency of the G — M20 time-scales. 

The close pair time-scales depend on the orbital decay 
times, with the smaller projected separations showing the 
greatest fractional variability. For 5 < Rproj < 20h~^ kpc, 
the observability time-scales are ~ 0.4 — 0.6 Gyr for the G- 
series simulations and ~ 0.1 — 0.3 Gyr for Sbc. These are 
slightly longer at 10 < Rproj < SOh~^ kpc, with time-scales 
~ 0.5 — 0.7 Gyr for the G-series sinmlations and ~ 0.2 — 0.5 
Gyr for Sbc. At larger projected radii < 50, 100 h^^ kpc, 
the highly radial orbit and the large pericentric distance 
parabolic orbit have significantly shorter (~ 0.3, 0.8 Gyr) 
and longer time-scales 1.4, 3.0 Gyr), respectively. The 
typical 10 < Rproj < 50h^^ kpc time-scales are ~ 0.7 — 1.2 
Gyr for the G3/G2 simulations and ~ 0.9 — 1.1 Gyr for 
the Sbc simulations. The typical 10 < Rproj < lOOft"^ 
kpc time-scales are ~ 1.9 Gyr for the G3 simulations and 
~ 1.1 — 1.4 Gyr for the Sbc simulations. Therefore the ob- 
servability time-scales for close pairs vary by a factor of 4-6 
at small projected separations and by a factor of 2-3 at larger 
separations. It is important to keep in mind that these are 
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Table 7. Equal-Mass Close Pair Timescales 



Simulation 


T(5 < Rproj < 20) 
(Gyr) 


T{10 < Rproj < 30) 
(Gyr) 


T(10 < Rproj < 50) 
(Gyr) 


T(10 < Rpro] < 100) 
(Gyr) 






Sbc-Sbc merj 


;ers 




SbcPPxlO 
SbcPPx4 


0.15 ± 0.18 
0.26 ± 0.21 


0.39 ± 0.22 
0.53 ± 0.31 


1.00 ± 0.17 
1.12 ± 0.17 


1.22 ± 0.20 
1.34 ± 0.19 


SbcPP 

SbcPR 

SbcRR 

SbcPPr- 

SbcPPr+ 

SbcPol 

SbcR 


0.15 ± 0.19 
0.27 ± 0.49 
0.16 ± 0.26 
0.13 ± 0.10 
0.26 ± 0.29 
0.10 ± 0.08 
0.08 ± 0.20 


0.35 ± 0.23 
0.54 ± 0.50 
0.44 ± 0.35 
0.43 ± 0.16 
0.50 ± 0.43 
0.34 ± 0.18 
0.18 ± 0.18 


0.90 ± 0.16 
1.08 ± 0.59 
0.97 ± 0.24 
0.78 ± 0.11 

1.38 ± 0.92 
0.92 ± 0.38 
0.36 ± 0.18 


1.20 ± 0.18 
1.37 ± 0.55 
1.31 ± 0.14 
0.99 ± 0.15 
3.06 ± 0.52 
1.45 ± 0.26 
0.84 ± 0.09 


SbcPPn=0 
SbcRn=0 


0.15 ± 0.20 
0.03 ± 0.06 


0.43 ± 0.28 
0.13 ± 0.19 


0.96 ± 0.17 
0.29 ± 0.21 


1.27 ± 0.15 
0.77 ± 0.05 


G-G mergers 


G3PP 
G2PP_^ 
GIPP^ 
GOPP^ 


0.39 ± 0.30 
0.43 ± 0.20 
0.58 ± 0.13 
0.67 ± 0.19 


0.72 ± 0.39 
0.60 ± 0.16 
0.52 ± 0.16 
0.39 ± 0.22 


1.21 ± 0.38 
0.71 ± 0.20 
0.64 ± 0.20 


1.85 ± 0.13 


G3PPn=0 
G2PPn=0'" 
GlPPn=0* 
GOPPn=0'= 


0.30 ± 0.38 
0.39 ± 0.16 
0.54 ± 0.09 
0.64 ± 0.17 


0.56 ± 0.39 
0.62 ± 0.17 
0.50 ± 0.17 
0.38 ± 0.21 


1.10 ± 0.37 
0.70 ± 0.20 
0.54 ± 0.20 


1.72 ± 0.13 



Rproj has units h ^ kpc. Timescales for simulations with starting separations less than maximum Rproj are 
not calculated. 



Initial separation is 70 /i ^ kpc. 
^ Initial separation is 56 kpc. 
Initial separation is 42 kpc. 



the observability time-scales for truly merging pairs {Tpair 
in Eqn. 16) and does not include the contamination correc- 
tion for non-merging pairs observed in projection ( p(merg) 
in Eqn. 16). Close pairs with large projected separations are 
more likely to be contaminated by non-merging galaxies, so 
the optimal separation distance is likely to be at intermedi- 
ate separations between 30-50 kpc. 

Unlike our simulated mergers where we know the 
merger parameters and initial galaxy properties a priori, 
it is generally impossible to recover these for each galaxy 
merger observed in large surveys of the distant universe. 
Ideally, one would like to determine an effective observ- 
ability time-scale for each method of identifying mergers 
in order to convert the number density of observed merg- 
ers observed into a galaxy merger rate. This effective time- 
scale should be weighted by the distribution of initial galaxy 
properties, mass ratios, and orbital parameters predicted for 
galaxy mergers by cosmological simulations. Our work here 
is a first step towards determining the mean observability 
time-scale and has concentrated on the systems most likely 
to be affected by dusty starbursts, i.e. gas-rich equal mass 
mergers of discs with small bulges. While this is significant 
improvement over previous estimates of the morphological 
disturbance time-scales, it is not sufficient to convert the 
observed fraction of morphologically disturbed and paired 



galaxies into a galaxy merger rate. Our next paper will ex- 
plore the merger observability time-scales of non-equal mass 
mergers needed to estimate the effective observability time- 
scale for a realistic population of mergers and calculate the 
galaxy merger rate. 

Galaxy mergers are often assumed to be associated 
with vigorous starbursts, such as observed for local ultra- 
luminous infrared galaxies (e.g. Sanders & Mirabel 1996). 
Our results here suggest that the timing of morphological 
disturbances can be offset from the peak in star-formation 
rate, especially if the stiff n = 2 supernova feedback pre- 
scription is correct. In general, the maximum morphological 
disturbances occur before the peaks in the star-formation 
rate. While asymmetries experience a sharp peak at the first 
pass lasting 100-200 Myr, the star-formation rate of the sys- 
tem remains enhanced above the initial rate for significantly 
longer after the first pass for both supernova feedback mod- 
els. Asymmetries also peak sharply at the final merger, while 
the star-formation peaks after the galaxies appear as single 
object when the nuclei coalesce. Therefore, the objects with 
highest star-formation rates may not always have the high- 
est asymmetries. This is particularly true for the stiff n = 2 
supernova feedback models, which experience a long burst 
at the final merger (~ 500 Myr) lasting longer than the high 
asymmetries (~ 200 Myr). During both first pass and final 
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Table 8. Equal-Mass Merger Remnant Properties 



Simulation 


G 


A/20 


C 


A 


SFiJ (Mo yr-l) 


Sbc-Sbc mergers 


SbcPPxlO (no dust) 


0.59 ± 0.03 


-1.75 ± 0.39 


4.1 ± 0.3 


0.05 ± 0.02 


5.82 


SbcPPxlO 


0.54 ± 0.01 


-1.80 ± 0.23 


3.4 ± 0.3 


0.15 ± 0.07 


5.82 


SbcPPx4 


0.55 ± 0.01 


-1.98 ± 0.07 


3.6 ± 0.1 


0.16 ± 0.05 


5.28 


SbcPP 


0.54 ± 0.01 


-1.93 ± 0.09 


3.2 ± 0.1 


0.00 ± 0.48 


5.24 


SbcPR 


0.53 ± 0.02 


-1.91 ± 0.11 


3.2 ± 0.1 


0.09 ± 0.02 


5.11 


SbcRR 


0.54 ± 0.02 


-1.77 ± 0.11 


3.1 ± 0.1 


0.14 ± 0.04 


6.78 


SbcPPr- 


0.53 ± 0.03 


-1.85 ± 0.20 


3.1 ± 0.2 


0.08 ± 0.10 


4.58 


SbcPPr-l- 


0.59 ± 0.04 


-1.94 ± 0.28 


3.5 ± 0.2 


0.01 ± 0.10 


3.68 


SbcPol 


0.55 ± 0.02 


-1.92 ± 0.14 


3.2 ± 0.2 


0.16 ± 0.05 


5.45 


SbcR 


0.52 ± 0.03 


-1.80 ± 0.22 


3.2 ± 0.2 


0.25 ± 0.04 


12.85 


SbcPPn=0 


0.59 ± 0.04 


-2.46 ± 0.38 


4.5 ± 0.6 


0.12 ± 0.08 


3.41 


SbcRn=0 


0.54 ± 0.04 


-2.21 ± 0.49 


3.9 ± 0.9 


0.22 ± 0.04 


8.33 


G-G mergers 


G3PP 


0.55 ± 0.05 


-1.93 ± 0.10 


3.2 ± 0.3 


0.03 ± 0.02 


3.11 


G2PP 


0.56 ± 0.02 


-1.83 ± 0.14 


3.2 ± 0.2 


0.01 ± 0.04 


1.97 


GIPP 


0.53 ± 0.01 


-1.92 ± 0.05 


3.2 ± 0.2 


-0.02 ± 0.03 


0.68 


GOPP 


0.50 ± 0.01 


-1.73 ± 0.06 


3.2 ± 0.8 


0.03 ± 0.04 


0.38 


G3PPn=0 


0.57 ± 0.02 


-2.37 ± 0.20 


4.2 ± 0.5 


0.04 ± 0.03 


1.88 


G2PPn=0 


0.57 ± 0.01 


-2.25 ± 0.12 


4.1 ± 0.4 


-0.01 ± 0.03 


0.95 


GlPPn=0 


0.60 ± 0.01 


-2.04 ± 0.07 


4.9 ± 0.7 


-0.05 ± 0.02 


0.13 


GOPPn=0 


0.52 ± 0.02 


-1.73 ± 0.06 


3.8 ± 0.5 


-0.01 ± 0.02 


0.03 



The properties of the simulated merger remnants observed 1 Gyr after the coalescence of the nuclei. 



merger, the majority of the enhanced star formation occurs 
in the nucleus which makes the quantitative morphologies 
appear more concentrated but not necessarily disturbed. If 
we ignore the effect of dust on the morphologies, G cor- 
relates directly with the star-formation rate. Because dust 
lowers the measured G value, the dusty Sbc final mergers 
do not show high G — M20 , while the less-obscured G-series 
final mergers do (Figs. [T2] [TT)) . 

Roughly 75% of the strongest starbursts in the local 
universe, ULIRGs, have quantitatively disturbed morpholo- 
gies, and two-thirds show multiple nuclei. Therefore the 
correlation between disturbed morphologies and peak star- 
formation rates appears to be better than what is implied 
by our models. We note that only one of our simulations 
(the gas-rich, radial orbit, isothermal feedback simulation 
SbcRn=0) reaches a star-formation rate at the final merger 
that is comparable to ULIRGs. ULIRGs may have higher gas 
fractions than our models, and often host active nuclei that 
could destroy or sweep out dust in the central regions during 
the final merger, both of which would tighten the correlation 
between disturbed morphologies and high star-formation 
rates. Finally, the star- formation rates of some ULIRGs may 
be overestimated, as some of the infrared luminosity may be 
from dust-heating by an AGN rather than star-formation. 
Any correlation between high star-formation rates and dis- 
turbed morphologies will depend on the star-formation indi- 
cator used to calculate the star-formation rate. Observed Ha 
luminosities and equivalent widths are expected to be high- 
est during the first pass and initial starbursts before the 
star-forming regions have been enshrouded in dust, while 



infrared luminosities will peak during the final merger after 
sufficient amounts of metals and dust have been produced 
(Jonsson et al. 2006). 

The vast majority of our equal-mass gas-rich merger 
remnants are decidedly disc-like and dusty, even those merg- 
ers which started with relatively low gas fractions. Previous 
studies of these and other equal-mass merger simulations 
have shown that the mass distribution of star particles follow 
r^''^ laws with steep central cusps consistent with or more 
concentrated than the light profiles of elliptical galaxies (e.g. 
Cox et al. 2006, 2008; Naab et al. 2006; Bournaud, Jog, & 
Combes 2005) and lie on the fundamental plane (Robert- 
son et al. 2006b). Only simulations with gas fractions > 
50% (Springel & Hernquist 2005; Robertson et al. 20Q6a) or 
mass ratios less than 1:3 (Naab et al. 2006, Bournard et al. 
2005) have been found to have merger remnants with mas- 
sive disc components. The masses of our merger remnants 
are dominated by the bulge component; our merger rem- 
nants appear disc-like and less concentrated than elliptical 
galaxies because we examine the g-band light profiles which 
are strongly affected by both dust extinction in the central 
regions and a bright but low-mass disc of young stars. Only 
models with less dusty remnants (as produced by n = 
supernova feedback) or models that ignore dust in the rem- 
nant produce remnants with high enough G, G, and low 
enough M20 values to be called spheroids. However, even 
these spheroidal remnants are forming stars at rates > 1 
M0 yr~^, and would not be classified as red E/SO or post- 
starburst E-l-A galaxies. If these simulations accurately rep- 
resent the end stages of the merger process, a number of 
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Figure 19. Left: G3 and Sbc remnants viewed from camO. Right: 
G3 and Sbc remnant viewed from cam4. All of the remnants show 
a low-mass dusty star-forming disc as well as a large bulge com- 
ponent. 



'green' Sb galaxies may be merger remnants (e.g. Hammer 
et al. 2005). However, our models may overestimate the ex- 
tinction and star-formation during the post-merger stages. 
Destroying dust by shock-heating would not be sufficient to 
produce true red and dead spheroids, and additional physics 
such as feedback from an active galactic nucleus may be 
needed to clear out the gas and kill star formation (see also 
Khalatyan et al. 2008). The remnants forming the fewest 
stars and with the least dust are those produced by the low- 
est mass systems. But the low-mass remnants are the least 
centrally concentrated of all the simulated remnants. Like 
their higher mass counterparts, the low-mass remnants have 
more gas (~ 10* M©) and star formation than typical dwarf 
ellipticals. 



6 SUMMARY 

We present a morphological analysis of a large suite of GAD- 
GET N-Body/SPH equal-mass gas-rich disc galaxy merger 
simulations which have been processed through the Monte- 
Carlo radiative transfer code SUNRISE. With the result- 
ing images, we have examined the dependence of quanti- 
tative morphology and projected separation in the SDSS 
g-band on merger stage, dust, viewing angle, merger ori- 
entation and orbital parameters, gas properties, supernova 
feedback prescription, and total mass. We have determined 
the time-scales for quantitative morphological disturbances 
in the Gini coefficient, M20, C, and A, and the time-scale 
during which close pairs lie at projected separations Rproj < 
20, 30, 50, and 100 kpc. We also examine the merger 
remnant morphologies and star-formation rates. 

• All of the equal-mass gas-rich merger simulations ex- 
perience quantitatively disturbed morphologies in G — AI20, 
G — A, and A at the first pass and/or the final merger. 
This is in good agreement with the morphologies and merger 
stages of the local ULIRG sample used to empirically cali- 
brate these quantities. However, merging galaxies observed 
between the first pass and final merger or after the coales- 
cence of their nuclei may not show disturbed G — AI20 and 
asymmetries. 

• The time-scale during which an equal-mass gas- 
rich merger may be identified is strongly dependent on 
the method used to find the merger. The G — M20 time- 
scales are the shortest of the morphological methods, but 
have the least dependence on the merger parameters with 
T{G - M20) ~ 0.2 - 0.6 Gyr. The asymmetry time-scales 
vary by a factor of 3-4 between 0.2 — 1.1 Gyr, and the 
G~ A time-scales are the longest, with T(G — A) 0.3 — 1.2 
Gyr. The close pair time-scales vary by factor of 2—6 with 
the orbital parameters, depending on the projected separa- 
tions adopted. At 5 < Rproj < 20h~^ kpc, the observability 
time-scales are ~ 0.1 — 0.6 Gyr. At 10 < Rproj < SOh~^ 
kpc, the observability time-scales are ~ 0.2 — 0.7 Gyr. At 
10 < Rproj < 50h~^ kpc, the typical observability time- 
scales are ~ 0.7-1.2 Gyr. At 10 < Rproj < lOOh'^ kpc, the 
typical observability time-scales are ~ 1.1 — 1.9 Gyr. 

• The presence of dust has strong impact on the quanti- 
tative morphological measurements, lowering G and G, rais- 
ing M20 throughout the merger, and raising A during the 
post-merger and remnant stages. 

• When dust is included, the time-scales for morpholog- 
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ical disturbances are most sensitive to the gas fraction of the 
merging galaxies, their pericentric distance, and relative ori- 
entation. The supernova feedback proscription and the total 
meiss of system do not significantly change the morphological 
time-scales. The relative orientations also affect the strength 
of the morphological disturbances, with prograde-retrograde 
and polar orientation showing the highest asymmetries. The 
timing of the disturbances also depends on orbital param- 
eters and gas fractions, with low gas fractions, large peri- 
centric distances, and highly radial orbits showing strong 
disturbances primarily during the final merger. 

• The timing of morphological disturbances is gen- 
erally offset from the peak in stax-formation rates, with 
strong morphological disturbances occurring before bursts 
of merger-induced star-formation and for shorter durations. 
Hence, not all merger-induced starbursts will exhibit mor- 
phological disturbances and vice versa. The mode of super- 
nova feedback and dust production also play important roles 
in the correlation between morphological disturbances and 
observed star-formation indicators. 

• The majority of simulated merger remnants observed 
> 1 Gyr after the coalescence of their nuclei appear disc- 
like and dusty in g-band light and are consistent with 
early-type spiral morphologies and star-formation rates. De- 
creased dust extinction would make most remnants appear 
more spheroidal, but would not affect the remnants' high 
star-formation rates (typically > 1 M0 yr~^). A major gas- 
rich merger without AGN feedback does not, by itself, pro- 
duce a red and dead spheroidal galaxy. 
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